Quantum kinetics and prethermalization of Hawking radiation 
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We reinvestigate the emission of Hawking radiation during gravitational collapse to a black hole. 
Both CGHS collapse of a shock wave in (l+l)-dimensional dilaton gravity and Schwarzschild collapse 
of a spherically symmetric thin shell in (3 + l)-dimensional gravity are considered. Studying the 
dynamics of in- vacuum polarization, we find that a multi-parametric family of out- vacua exists. 
Initial conditions for the collapse lead dynamically to different vacua from this family as the final 
state. Therefore, the form of the out- vacuum encodes memory about the initial quantum state 
of the system. While most out- vacua feature a non-thermal Hawking flux and are expected to 
decay quickly, there also exists a thermal vacuum state. Collectively, these observations suggest an 
interesting possible resolution of the information loss paradox. 

PACS numbers: 



I. INTRODUCTION 

Hawking' s discovery of black hole radiance [l[ in 1974, 
along with Bekenstein's conjecture [2| of the generalized 
second law of thermodynamics in 1973 inspired a con- 
ceptually pleasing connection between general relativity, 
thermodynamics and quantum theory. It also led to a 
serious conundrum. As Hawking showed, when a quan- 
tum field theory is put on a fixed black hole background, 
an observer at spatial infinity will detect a flux of the 
quanta of that field in a mixed quantum state described 
by a thermal density matrix. This thermal flux leads 
in time to the complete evaporation of the black hole. 
Thus, if matter quanta are thrown into an already- formed 
black hole, any quantum information encoded in their 
wave function is lost to the outside once those quanta 
cross the horizon. Because the Hawking radiation itself 
encodes no information other than the size of the black 
hole, the in-fallen information cannot be recovered as the 
black hole evaporates. Thus, once the black hole evap- 
orated completely, a pure incoming quantum state has 
been transformed into a mixed outgoing quantum state 
in apparent contradiction with the unitarity of quantum 
mechanics. This is the essence of the famous information 
loss paradox 0, Q . 

The black-hole information loss paradox is one of the 
most pressing problems at the interface between quan- 
tum mechanics and General Relativity, and it does not 
seem that one can appeal to quantum gravity for a solu- 
tion. The quanta of the Hawking radiation are produced 
in the vicinity of horizon, where the gravitational field 
is weak for macroscopic black holes; therefore, the ap- 
proximation of doing quantum field theory on a fixed 
curved background spacetime should be valid. Neither 
can one expect the paradox to somehow be resolved once 
the evaporating black hole becomes of Planckian size and 
quantum gravity does kick in. There is far too much in- 
formation encoded by the in-falling matter of a stellar 



black hole to be stored in a Planckian remnant - the in- 
formation must be released long before quantum gravity 
should matter. 

In an attempt to understand the information loss prob- 
lem in the context of quantum gravity, Callan, Gid- 
dings, Harvey and Strominger (CGHS) studied (1 + 1)- 
dimensional dilatonic black holes The major ad- 

vantages to the CGHS model are that, unlike the 4- 
dimensional Einstein gravity case, (a) (l + l)-dimensional 
quantum gravity is a renormalizable theory and (b) the 
expectation value of the energy-momentum tensor of the 
Hawking radiation is directly related to the conformal 
anomaly of the theory [6] and is therefore explicitly cal- 
culable. This allows back-reaction to be included into 
the model relatively easily. Once this is done, informa- 
tion still appears to be lost |7j due to the emission of 
thermal Hawking radiation. 1 

There have been several classes of proposed resolutions 
of the information loss problem. Most (if not all) of these 
proposals [|| seem to bring with them other physical or 
conceptual difficulties: 

(1) Information really is lost [3|, implying that our 
present laws of quantum mechanics are insufficient to 
address the formation and evaporation of a black hole. 
Thus one either has to abandon the idea of energy con- 
servation by introducing a non-unitary time evolution or 
allow quantum information to dissipate and therefore ul- 
timately violate unitarity of quantum mechanics (see for 
example P3). 

(2) Information is not lost; instead, it is stored inside 
the black hole. The black hole does not evaporate com- 
pletely, but leaves behind a remnant that has a mass on 
the order of Mp [ll| . This remnant contains all the origi- 
nal information. There are two major problems with this 
idea: (a) it requires infinitely many species of remnants 
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1 Recently, evaporation of 2-dimensional CGHS black holes with 
backreaction effects taken into account has been revisited in 
revealing a non-thermal power-law-like spectrum of Hawking ra- 
diation. 
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(since there are infinitely many initial states leading to 
gravitational collapse), which, in addition to aesthetic 
objections, leads to an unacceptably high rate of pair 
production of such remnants in a gravitational field [12j; 
and (b) there seems to be no mechanism to stabilize such 
remnants. A related idea is that a "baby universe" is 
formed within the black hole and the in-falling informa- 
tion gets transferred to that universe. 

(3) The physics of a black hole and of Hawking radia- 
tion is non-local (see for example [l3|, [HI)- Deviations 
from locality would automatically imply corrections to 
the Hawking result. One way to introduce this non- 
locality is to assume that the degrees of freedom in- 
side a black hole are duplicated in the exterior [HI, [HI . 
The information thereby remains intact in the exterior. 
This however seems to be at odds with the superposi- 
tion principle of quantum mechanics [l7j. Still, the idea 
of non-local physics somehow resolving the information- 
loss paradox presently seems to be the most popular in 
the literature. 

(4) Information escapes during collapse and evaporation. 
This can happen if information is subtly encoded into 
correlations among the particles that fall into and escape 
the black hole or if, contrary to Hawking's calculation, 
the Hawking radiation is not completely thermal [l8| and 
can carry out small amounts of information. It is pos- 
sible that (4) is automatically included in (3), once we 
learn how to describe the non-local physics of gravita- 
tional collapse. 

In this paper we try to address both possibilities (3) 
and (4) for the scenarios of CGHS collapse as well as 
(3 + l)-dimensional spherically symmetric Schwarzschild 
collapse according to the following program: 

• All the calculations will be done in the 
Schwarzschild slicing, where the collapse takes 
an infinite amount of time to complete but the 
collapsing shell of matter approaches its apparent 
horizon exponentially quickly. 

• Following Hawking, we shall quantize the scalar 
field f(t,r) on the fixed background spacetime ap- 
propriate to the collapse, neglecting back-reaction 
effects. The difference from the standard ap- 
proach is that we shall use Heisenberg quantiza- 
tion, where the modes of the field are considered 
time-independent while their time-dependent am- 
plitudes are promoted to quantum-mechanical op- 
erators. 

• Using the results of Heisenberg quantization, we 
shall study the dynamics of in- vacuum polarization 
in order to find how and when the out-vacuum of 
the scalar field is formed. We shall also determine 
the explicit form of the out- vacuum. 

• Finally, we shall calculate the Bogolyubov coeffi- 
cients between the in- vacuum and the formed out- 
vacuum in order to find the spectrum of Hawk- 
ing radiation. We shall also calculate the response 



function of a detector carried by a fiducial observer 
located at spatial infinity. 

As we shall see, this program leads to several surpris- 
ing results. We find that an out- vacuum is formed as 
a Bose-Einstein condensate of "in" excitations within a 
time interval St ~ Rs after the onset of the collapse. In- 
stead of a unique out- vacuum, there exists a continuous 
family of "out" vacua. Which particular vacuum from 
this family becomes the post-collapse final state depends 
on the initial state of the collapse process. A calcula- 
tion of the response function of a detector shows that, 
for most of these vacua, the Hawking flux is not thermal 
(and divergent), while for one particular vacuum from the 
family it is thermal (and finite). This fact has important 
implications for the information loss paradox. 

The paper is organized as follows. In Section [H] we 
briefly review the classical gravitational collapse of a thin 
shell of matter in both (l + l)-dimensional dilaton gravity 
(Section III Aft and (3 + l)-dimensional general relativity 
(Section III C|) . We perform Heisenberg quantization of 
the scalar field on a fixed collapsing CGHS background 
in Section IIIII To check that the quantization was per- 
formed correctly, we assume a particular form for the out- 
vacuum and calculate Bogolyubov coefficients between 
in- and out- vacua in Section [IVJ The calculated Hawking 
flux reproduces the classical result by Giddings and Nel- 
son [l9j. Then, in Section [V] we derive the quantum ki- 
netic equation describing Bose-Einstein condensation of 
"in" excitations and the process of formation of the out- 
vacuum. We find the explicit form of the out-vacuum. 
The same analysis is repeated for (3 + l)-dimensional 
spherically symmetric Schwarzschild collapse in Section 
IVII We summarize our results and muse extensively on 
the resolution of the information loss paradox in Sec- 
tion IVIII In particular, subsection IVII Al contains the 
review of our major results. Physical content of quan- 
tum kinetic equations discussing particle production in 
QFTs on fixed curved background is discussed in subsec- 
tion [VIIBl Properties of out- vacua and physical observ- 
ables are discussed in subsection IVII CI and implications 
of our results for information loss paradox — in sub- 
section [VHDj Finally, the collapse of a thin spherically 
symmetric shell on the formed primordial black hole is 
considered in Appendix O while the case of two con- 
centric collapsing shells is discussed in Appendix [Bl In 
Appendix [C] we explain how our results help to resolve 
the famous issue of WKB breakdown in the near horizon 
regime for solutions of the Wheeler-de Witt equation de- 
scribing collapse and Hawking evaporation. 

II. GRAVITATIONAL COLLAPSE IN 
CLASSICAL THEORY 

We first review some of the aspects of the classical so- 
lution describing a thin shell collapsing to a black hole in 
two-dimensional dilaton gravity [5| and four-dimensional 
Einstein gravity (with spherical symmetry). Our major 
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concern will be to compute the dynamics of collapse in 
Schwarzschild-like coordinate systems where the appar- 
ent horizon, as seen by a fiducial observer at spatial in- 
finity, is spacelike, full collapse time (the time before the 
collapsing body /shell reaches its event horizon) is infi- 
nite, but the approach of the shell to its apparent horizon 
is exponentially rapid. 



A. CGHS collapse: classical solution 

Two-dimensional dilaton gravity is the theory de- 
scribed by the action 



d 2 x 



^[e-^(i? + 4(V(/)) 2 +4A 2 )-(V/) 2 ], (1) 



where g, </>, and / are the metric, dilaton, and matter 
fields, respectively, while A 2 is a cosmological constant 
(see |20j for comprehensive review). The equation of mo- 
tion for the dilaton <\> is given by 

= 2e- 2 * (V„V^ + (VV + A 2 - (V0) 2 )) , (2) 

while the dynamics of spacetime is determined by 

= 2e~ 2 ^ (4V 2 + 4A 2 + R - 4(V0) 2 ) . (3) 

It is convenient to analyze the theory in conformal 
gauge and introduce light-cone coordinates x ± = x° ±x 1 . 
Then, the metric takes the form 



ds 2 = —e 2p dx + dx 
with metric coefficients 

1 



2p 



9+- = 9-+ 

The Ricci scalar R present in (j3J) is given by 
R = 8e- 2p d+d-p. 



(4) 



(5) 



(6) 



In terms of p, <j> and /, the equations of motion can 
then be written as 

= e - 2 (*+rt (A 2 e 2 " - 4d+d-(f> + 2d+d- P + 

+ Ad+<j>d-<t>) (7) 

=e- 2 ^ (2d+d-<p - 4d+4>d-4> - \ 2 e 2p ) (8) 

=d + d.f, (9) 

The c/++ and g components give constraint equa- 
tions 



-20 



{4d±pd± 



2d 2 ±(j 



'^(d±f) 2 - (10) 



Since the gauge is not fixed completely by the condition 
(|5j), we can choose the Kruskal condition 



<t> = P, 



(11) 



so that the metric (j4]) reduces to 

ds 2 = -e 2 ^dx + dx~. 

The Eq. (j8]) now acquires the form 

d+d-(e- 2 *) = -A 2 

and upon solving yields 

e" 2 " = ^ - \ 2 x + x~ 
A 



(12) 



(13) 



(14) 



where M/X is an integration constant. (The linear terms 
in x ± are omitted, since they represent only constant 
shifts in those coordinates.) M is interpreted to be the 
mass of an eternal black hole. In the limit M = 0, this 
solution is known as the linear dilaton vacuum. 

If one starts with the linear dilaton vacuum and then 
sends a pulse of the matter scalar field /, a black hole is 
produced. For a sharp left-moving pulse, we can write 



1 



s(x + - 4). 



(15) 



From ([TO]) , again choosing p = </>, the general solution for 



-2cf) 



is given by 



" 2 * = -X 2 x+x 



-X 2 x + , 



J dx ,+ ^ dx"+5(i 



."+ ^+ 



Xx, 
M 
Xx£ 



4) 



(x + -x+)G(x+ -£+), (16) 



where Q(x — y) is the Heaviside theta function. Thus for 
x + < Xq this solution corresponds to the linear dilaton 
vacuum solution, while for x + > Xq it describes the black 
hole of the mass M . The spacetime metric takes the form 

dx + dx~ 



ds z 



-X 2 x+x~ - t^(x+ - x£)G(x+ - x+) 



(17) 



B. CGHS collapse: coordinate systems 

Here we will discuss the solution describing the (1 + 
l)-dimensional collapsing shell in three different co- 
ordinate systems: Kruskal, Eddington-Finkelstein and 
Schwarzschild (2lj . 

According to the expression ([TT|h the solution has dif- 
ferent forms in the regions x + < Xq and x + > Xq . 



1. x + < Xq region 

Inside the collapsing shell one has e~ 2( ^ 
so that the metric is 



-X 2 x+ r j 



ds 2 = 



dx + dx 
X 2 x+x~ 



(18) 



This coordinate system is a direct analogue of the 
Kruskal-Szekeres coordinate system for the (3 + 1)-D col- 
lapsing spherically symmetric shell. For a more physical 
interpretation, one can go to the Eddington-Finkelstein 
coordinates, which are defined by x ± = ±e ±Xa . The 
metric then is given by 



ds 2 = —dander , 



(19) 



so that the spacetime inside the collapsing shell is just flat 
Minkowski spacetime. For the Eddington-Finkelstein co- 
ordinate system, the inside region corresponds to Aa + < 
ln(Xx+). " 



Singularity 




FIG. 1: 
setup. 



Penrose diagram describing the CGHS collapsing 



Writing cr + = t a + r a , we find that the region x + < Xq 
corresponds to r a < A _1 ln(Axo") — t a . The equation 
r a = A -1 \u(\xq) — t a gives the position of the moving 
shell from the point of view of an observer inside the 
shell. According to such an observer, the shell moves at 
the speed of light. 



The region x + > Xq in the Kruskal-Szekeres coordinates 
corresponds to A<r + > ln(AxJ) or r a > A -1 Ih(Xxq) — ta- 
in the Eddington-Finkelstein coordinates. 
The scalar curvature is given by 



R = 



4MA 



4MA 



M/X - X 2 x+x~ M/X + e 2Ar - ' 



(22) 



where r a = (<j + — a )/2. 

As we see, a singularity is reached at x + x~ 



t 2 



x 2 = in the Kruskal-Szekeres coordinates or at r a = 
log(— M/X) in the Eddington-Finkelstein coordinates. 
One would naturally expect this singularity to be hid- 
den by the event horizon and this is indeed what hap- 
pens in this setup [I?}- If we assume for a moment 
that we are dealing with a (3 + l)-dimensional spherically 
symmetric setup instead of (1 + l)-dimensional one, the 
area of 2-spheres (corresponding to points in the (1 + 1)- 



dimensional setup) would be given by % e 



-2cf)(x + ,x~) 



Based on the form of the solution for the dilaton, it is 
easy to check that all the 2-spheres (points) such that 
x~ < —M/X 3 Xq = xh, x + < Xq are trapped. The ap- 
parent horizon is the boundary of such region given by 
x~ = xh- For the given setup, it coincides with the 
event horizon, so we shall denote it further simply as the 
horizon. 

In the Eddington-Finkelstein coordinates, the horizon 
is located at <r + = — oo (or a~ = +oo). From above, it 
is only reached as 



-A(2£ CT -A -1 ln(Ax+)) 



Xx^e 



-2\t a 



o, 



(23) 



that is, at t a — >• oo. 

As we see, the BH formation takes infinite time in the 
Eddington-Finkelstein coordinates, but the approach to 
it is exponentially rapid. Since we would like to make 
parallels with a spherically symmetric (3+l)-dimensional 
gravitational collapse, this class of coordinate systems is 
naturally the only one that interests us. 

Another important coordinate system where the 
dynamics of collapse has similar features is the 
Schwarzschild-like one. Let us define 



A 



2. x + > Xq region 



In the region outside the shell one has e 

2 ^y,+ n 



-2cf) 



M/X 



in the Kruskal-Szekeres coordinate system, so 



that the metric is 



ds 1 



dx+da 



M/X - X 2 x+x~ 



(20) 



The Eddington-Finkelstein coordinates are now defined 



by x^ 



and e 



-X(x~ + M/X 3 ), and the 



metric is given as 



ds 2 



dander 



1 + (M/A)e A (° 



(21) 



1 (M 
2A ln U 



„2Ar c 



1 fe\ 
2A ln U 



so that the metric can be rewritten in the form: 

dr 2 



(1 - e 1 " 2 ^) dtl 



1 - e l"2Ar- 



(24) 



(25) 



As usual, there is a coordinate singularity at r = 
(2 A) -1 = Rs, corresponding to the location of the hori- 
zon. 

Let us determine the equation for the motion of the 
shell in this coordinate system. Recall that the equa- 
tion r a = ^-log(Axo") — t a describes the motion of 
the shell in Eddington-Finkelstein coordinates; hence in 
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Schwarzchild-like coordinates, we find that 
R(t) = R S + Rs log (l + ( e (Ro-Rs)/Rs _ i) e -ta/Rs^ 
^R s (l + ( e (Ro-Rs)/Rs _ i)e-WHs) (26) 



where the position Ro of the shell at t = is given by 
e (R -R s )/Rs _ i — ( Xx o ) anc [ i n ^ e second line we have 



used the late-time (near horizon) limit t G ^> Rs. Again, 
we can see that the shell reaches the horizon only as 
t a — > oo. 

Finally, at early times t a — )• — oo the spacetime is flat 
up to exponentially small corrections. Indeed, from the 
Eq. ((26} we find that 



R(t) ~ Rq — t a 



(27) 



assuming that Ro ^> Rs (this apparently holds for suf- 
ficiently large shells). The shell moves with the speed 
of light, the spacetime inside the shell (at r < R(t)) is 
Minkowski while the metric of spacetime outside it (at 
r > R(t)) is given by the Eq. (|2S)). Since R(t) > R > 
Rs at early times, we come to the conclusion presented 
above. 

The overall physical picture of collapse described above 
can be fully represented by the Penrose diagram on the 
Fig. [D depicting the causal structure of the collapsing 
classical (1 + l)-dimensional spacetime. 



C. (3 + l)-dimensional Schwarzschild collapse: 
classical solution 

The case of (3 + l)-dimensional spherically symmet- 
ric classical gravitational collapse of an infinitely thin 
shell is in many respects similar to the case of the two- 
dimensional CGHS collapse. This should not be of any 
surprise: if the metric A gik of the 4-dimensional collaps- 
ing spacetime is represented in the form 



-2cf) 



9ik = 9ik + 



(28) 



where gik is the metric on the r — t plane, Sik is a unit 2- 
sphere metric, k is a constant with dimension of inverse 
length, and angular degrees of freedom are integrated 
out, the Einstein-Hilbert action of the 4-dimensional 
problem reduces to the one of two-dimensional dilaton 
gravity (see, for example, [8[). 2 



2 Strictly speaking, simple equivalence between 2- and 4- 
dimensional spherically symmetric collapse problems is of course 
lost at the quantum level. In particular, the form of the one-loop 
effective action responsible for the trace anomaly is different for 
the two cases p2l. [23j. Nevertheless, the results which we present 
below are equally valid for both 2-dimensional and 4-dimensional 
spherically symmetric collapse problems as should be clear from 
the procedure of their derivation. 
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FIG. 2: Penrose diagram describing 4-dimensional spherically 
symmetric shell collapsing into Schwarzschild black hole. 



The spacetime inside the collapsing spherically sym- 
metric shell of radius R = R(t) is Minkowski due to 
Birkhoff 's theorem: 



ds 2 = -dT 2 + dr 2 +r 2 dn 2 , 



(29) 



while outside the shell (r > R(t)) it is described by the 
Schwarzschild metric: 



ds 2 



Rs 



dt 2 



r 



dr 2 +r 2 dn 2 . (30) 



The exterior metric has a singularity at r = covered 
by the spacelike horizon at r = Rs, but the latter is 
not formed as long as the position of the shell R(t) > 
Rs. It is well known that from the point of view of a 
fiducial observer at spatial infinity the (3+l)-dimensional 
gravitational collapse requires an infinite amount of time 

ism. 

Indeed, let us derive indeed the equations of motion 
for the shell. To properly match the timelike coordinates 
inside and outside the shell, one has to refer to the proper 
time of an observer living on the shell. One finds [26| , [27j : 



(31) 



where B(t) = 1 - ^fy . 

The equations of motion for the shell are derived from 
the effective action [23 




eff 



4tto- / dTR 2 



dRY 
dTJ 



2irGaR j , 

(32) 
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where a is the mass density of the shell. One finds 
B 3 / 2 R 2 H 



B 2 ~ (f ) 



2 4trt(1 - 2irGaR s ) ' 



(33) 



where H is the conserved (total) mass of the shell. In the 
near horizon approximation the solution of this equation 
reads 



R(t) ~Rs + (Ro ~ Rs) exp (-t/Rs). 



(34) 



The expression ([34]) is close to ([26]) , the one describing 
the case of a (1 + l)-dimensional shell collapsing into 
CGHS black hole. 

It is clear what happens physically at early times t — »• 
— oo. The shell moves with the speed of light, and as 
long as R(t) ^> Ro ^> Rs, the spacetime is Minkowski 
up to poly normally small corrections. 

Once again, the classical causal structure of the colaps- 
ing spacetime can be represented by the Penrose diagram 
on the Fig. [2] featuring the event horizon. 



III. CGHS COLLAPSE: MODES, THE ACTION 
AND ITS DIAGONALIZATION 



After the brief review of (1 + 1)-D and (3+l)-D classical 
collapse, we are now in the position to discuss quantum 
effects and their influence on the collapse scenario. 

In this Section we consider quantization of a massless 
scalar field / on the fixed collapsing CGHS background 
([25]) . Equation of motion for the field / reads as usual 



V 2 /: 



1 d 



g dx % 



99 



df_ 

dx k 



0, 



(35) 



where is the metric ([25]) . The general solution of ([35]) 
can be represented as a superposition of eigenmodes of 
the V 2 operator: 



/ 



f(t,r)= / dk-a k (t)f k (r) 



(36) 



To quantize the field /, we then promote the time- 
dependent amplitudes ak(t) into quantum mechanical op- 
erators. Note that such a prescription is a bit unusual. 
Instead of the usual procedure of the Schrodinger quan- 
tization of the field f(t, r), i.e., representing the field as 
a superposition of time- dependent modes 



f(t,r) = dk- a,kfk(t,r) + 



c.c. 



(37) 



with a/c, a\ being time-independent canonical creation 
and annihilation operators, we cast the time depen- 
dence on the amplitudes a(t) themselves. The main 
reason for our choice is convenience — since the col- 
lapse problem is intrinsically time-dependent (the posi- 
tion R(t) of the shell is time-dependent and only reaches 



the Schwarzschild radius at t — » oo in the slicing which 
we use), it makes sense to use a formalism where the 
consequences of this time dependence can be found more 
easily. The major drawback of our quantization formal- 
ism is that the relativistic invariance of observables is less 
transparent. 

To determine the exact form of the eigenmodes /fc(r) 
we observe that (a) the modes are naturally divided into 
left-moving and right moving ones and (b) the form of 
the modes /fc(r) is different in the regions r > R{i) 
and r < R(t), where R(t) is the position of the shell. 
Naturally, the questions of normalization and matching 
the modes between interior and exterior regions become 
important. First, it is impossible to match the time- 
independent modes inside and outside the shell for an 
arbitrary moment of time t and r = R(t), since one 
would have to match the coefficients dk(t) at arbitrary 
t as well. However, after the quantization is performed, 
the coefficients a^it) carry the information about the oc- 
cupation numbers for a given mode /fe(r), 3 and these oc- 
cupation numbers do change with time (particles are pro- 
duced by the changed background spacetime). Instead, 
we shall match the modes only at the initial moment of 
time to assuming the vacuum-like initial state for them 
at to — » — oo. It is also natural to choose the space- 
like hypersurface R > Ro for the normalization of modes 
(note that if Ro ^> Rs, the spacetime is nearly flat — 
it is flat inside the shell, while all curvature invariants 
are very small outside the shell, so one would just have 
a usual Minkowski normalization for the modes). 

Taken these considerations into account, we find for 
a given mode fk with momentum k (the index L corre- 
sponds to left-moving modes, while the index R — to 
right-moving ones) 



fk,L/R(r) 



1 



(38) 



in the exterior of the shell, i.e., at r > R(t), where R(t) 
is the position of the shell ([26]) and 



fk,L/n(r) 



2k 



(39) 



inside the shell. 



To quantize the matter field /, we shall first need the 
relationship between the internal time, denoted here as 
T and the external time, denoted t. These can be found 
from the proper time of an observer "riding" the shell, 



3 For example, in the WKB regime, when occupation numbers are 
large, one can roughly estimate a&(t) ~ sj^k times the time- 
dependent quantum phase. 
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denoted r: 



dT 



dt 



-if 



so that 



where from (|2 



F 



F = 1 — e l ~ R / R s 



(40) 



(41) 



(42) 



To find the behavior of T as Rq- — »• ii^s we need to know 
i?. It is easy to show that the solution for in ([26]) 
results exactly to 



(43) 



Correspondingly, substituting (|^3|) into flH| gives T = 
F as i? -> 

Now, following [27], the action of the scalar field / can 
be written as the sum of two parts corresponding to the 
regions inside and outside the shell. One has 



S — Sin H~ S 



where 

Sin : 

and 

S ut 



J dt i ■ 



I di f 

J JR{t) 



dr 



dr 



F 



F(d r f) : 



(44) 



(45) 



(dtf) 2 



1 _ e l-r/Rs 



"(1 



Kdrff 



(46) 



Let us find out what kind of approximations simplifying 
the action ([44]) are possible. 

Since only the time dependent amplitudes a& are con- 
sidered quantum mechanical variables, after substituting 
the expansion ([36]) into the action we can explicitly take 
integrals over r. We then find 



S 



dtdkdk' 



-Kkk' — -Gkk' 



(47) 



where the kinetic energy contribution is 

Kkk' =dk,hM^h aa dk>L + dk,RMu3 )aa ak'R+ 



■a kiR M kk r a k ,, L + a } kL M kk / a k >, L - 



• t Ti/rRR,ca- i -t -\/rRL.ca- , 

a k,R M kW a k>,R + ak,R M kk> a k>,L + 

c.c., (48) 



while the gradient energy contribution is 

Gi\rLL,CLCL , -\rRRyCLCL • 

kk' =a>k,L±Vkk' a k' ,L -T Uk,R^kk' a k' ,R + 

■ t\tRL ,aa i t -\rLL.ca , 

+ a k,RN kk r aki,L + a kL N kk ,> a k ,, L + 

, t AT RR,ca , t AT RL,ca , 

+ a JU J W a k',R + a k,R N kk> a k >, L + 

-h c.c. (49) 

The kernels M^/ and 7V^/ are defined according to 



j\jLL,aa 

kk ' F(t) 

r+oc 
'R(t 



l r 

Jt) Jo dr f^ L ^f k ^ L ^ r ^ 
f+°° dr 

/ 1 iZ^fk,L(r)fv,L(r), (50) 



1 r R W 

jyjRR^aa ^ I 



drfkM r )fk'Mr) + 



/ 1 1 - r/R J kA r )fk'A r )^ ( 51 ) 



,l*(t) 

M ™=Tv7V / dr/ fc|L (r)/ fc .,fl(r)+ 



Jo 



etc., 



^"• aa = / dr (l - e 1 -^) fU(r)f'k>,L( r )+ 

JR(t) V 7 

+F(t) [ m drf k , L (r)f' k , tL (r), (53) 

^• aa = f° dr (l - e W* S ) /^(r)/^, fl (r)+ 

+F(t) / drf KR {r)f k ,^{r), (54) 

JO 

A^ aa = f° dr (i - e wM 

«/ R{t) 

+F(t) f m drfi L (r)f' k , }R (r), (55) 

«/ 

etc. 

Let us explicitly calculate (LF, aa^) and (i?i?, aa^) ma- 
trices. The matrix M kk ,' aa has two contributions - 
from the regions inside and outside the shell. The first, 
interior one, is 

I rR(t) 



F(t) 



[ drf k ,L(r)fk>A r ) 
Jo 



R(t) 



2F(t)Vkk' Jo 
H(k-k',R{t)) 



dre i{k-k') ^ 56 ) 
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where we have introduced the function 



H(k-k',x) 



bo ( ( 



z(jfe-jfe') 



t i(k—k f )x 



(57) 



The second, exterior, contribution to M^ ,aa can be in 
turn expressed as 

dr 



R{t) 



I _ e l-r/R s 
1 



A,L(r)/^ |L (r) = -*(*-*/)- 
-.H(k - k\ Rs log (e m/Rs - e)) (58) 



- L S (k - k') - H(k-V,(*x+)*e-*/ R s) 
2k [ } 2Vk& 

Using the properties of the function H(k — one can 

further simplify these expressions. As k — k' — >• (or x —> 
0), the function if(fc — fc^x) approaches the finite value 
1. At sufficiently large k — k' and x it rapidly oscillates 
effectively cutting off the correspondent contributions to 
the action at \k — k'\ ~ \jx. Therefore, integrals of the 
form j dkdk'H(k — k 1 \x)A k A^ can be estimated as 



J dkdk f H{k-k f ,x)A k A%, 



dkdAkH(Ak, x) A k A* k _ Ak 



(59) 



A k \\ 



where Ak = k — k' and we used the quasiclassic ap- 
proximation assuming that the amplitude A k is a slowly 
changing function of the momentum k. 

The Eq. (|59|) implies that the function H(k — k f : x) 
can be safely substituted by the delta function S((k — 
k')x) = ^S(k — fc'), when x is sufficiently large. How 
large is "sufficiently large"? As we shall see soon, al- 
though the "in" quantum state of the field f(t,r) (at 
t — >• — oo ) is described using momentum quantum num- 
bers fc, the relevant quantum numbers for the "out" state 
(at t — )• +oo) are not fc, but exponentially redshifted mo- 
menta p ~ k- exp (t/2Rs). Therefore, while the function 
H(k — k',x) does not completely coincide with ^S(k — k') 
being a smooth non-singular distribution smeared over 
the interval \k — k'\ ~ 1/x of "in" momenta, its width in 
terms of "out" momenta \p — p'\ ~ exp (t/2Rs)/x van- 
ishes as t — )• + oo. 

Keeping this in mind, we find that the interior con- 



fer \k - k'\ <C l/R(t) ~ l/Rs at late times t -> +oo. 
The interior contribution is approximately diagonal for 
all k and k' at early times. On the other hand, it is 
only diagonal for k^k' ^> l/Rs, while all the modes with 
fc, k' < l/Rs mix between each other. 4 We also note 
that the prefactor 1/F{t) ~ 1 at early times t — >> — oo 
and blows up exponentially at late times t +oo. 
Similarly, the exterior contribution to the matrix 

Mk£ ,aa is found to be 



R(t) 



dv 

^^-r/Rs f^m.dr) « Yk Kk - k>) (62) 



in the asymptotic future t — >• + oo, while in the asymp- 
totic past this contribution vanishes entirely. Taking into 
account the behavior of F(t), we therefore conclude that 



LL,* 
kk' 



1 

2k 



S(k - k') 



at t — >• — oo, while at £ — >• +oo it is 



LL,i 
kk' 



2kF(t) 



S(k - k') 



(63) 



(64) 



for the modes with momenta fc, k' ^> l/Rs- In principle, 
these are the modes of the main interest for us as the 
geometric optics limit can only be applicable for such 
modes. 

A similar analysis of the terms contributing to the ma- 



trix N, 



shows that 



LL,aa 
kk' 



k 



S(k - k') 



(65) 



for both asymptotic past and future infinities t — >• d=oo. 

Corrections to the expressions (|63|) . ([64]) and ([65]) are 
suppressed by the parameter exp (—\t\/Rs), very small at 
t — >• ±oo. At intermediate times t — >• these corrections 
become of the order one, but repeating our analysis and 
using the properties of the function H(k — k'^x) suggests 
that the only modes affected by the change of background 
spacetime between t — » — oo and t ~ are deep IR modes 
with k,k' < 1/Rq, where R = R(t = 0) > R S - 5 If we 
are not interested in those much, we can just set the 
initial quantum state ^(a kl t = 0) of the system at t = 
as the direct product of the "in" (Minkowski) vacuum 
states for the modes with different k. 

Finally, it is worth noting again that at t — >• +oo the 



matrix M, 



LL,aa 
kk' 



is not diagonal for the modes with mo- 



tribution to the matrix M^' aa is approximately given menta M' < V^s- Therefore, these IR modes with 



by 



1 



F(t) Jo 



R(t) 



drf k , L (r)f k , yL (r) 



2kF(t) 



S(k - k f ) (60) 



at early times t — >• — oo, when the size of the shell R(t) 
is very large, and 



1 



rR(t) 



F(t) Jo 



drfkAr)fk' L (r) 




(61) 



4 In fact, the modes with arbitrary /c, k' mix between each other 
within the interval \k — k'\ ~ l/Rs but for k, k' ^> l/Rs this 
interval is small compared to the values of k and k' . 

5 As will become clear from the further discussion, this approxima- 
tion is rather rough but sufficient to address our goals, namely, 
to determine the structure of the out-vacuum and calculate Bo- 
golyubov coefficients between "in" and "out" vacua, see Section 

[VTT1 
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fixed momenta do not diagonalize the action ([44]) and 
the associated Hamiltonian. As t — >> +00, these modes 
redshift away (as we shall see soon, the physical frequency 
at t —> +00 is k^/F(t)) and do not affect the observables. 

Applying similar considerations to other terms con- 
tributing to the action of the scalar field / and focusing 
only on the modes with k ^> 1/Rs, we conclude that the 
action for a given mode can be written approximately as 



S 



dk ■ k\dk\ 



(66) 



where F(t) is a smooth function with the asymptotic be- 
havior F(t) —> 1 at t — » —00 and F(t) — » at t — » +00. 
In particular, choosing the function 



F(t) 



1 

-t/R s 



when t < 
when t > 



(67) 



is sufficient to capture all the essential physics of the 
problem — as should be clear from the previous dis- 
cussion, changing F(t) while keeping its asymptotic be- 
haviour intact only affects the physics of modes with 
k < 1/Rq < 1/R S . 

The field ak entering the action ([66]) is not canonically 
normalized. In order to canonically normalize it, one has 
to introduce the change of variables 



A k = 



Vk' 



(68) 



Correspondingly, the modes ([38]) . ([38]) have to be changed 
to 



fk(r) = Vkf k (r) 
and acquire the form 

1 n ±ikr 

(e 2Xr - e) 



(69) 



Mr) 



-t=e ±lkr if r < R s 



V2 



In what follows we always consider the modes fk and 
operators A^, not the modes fk and operators a& and 
shall therefore simply omit the tilde. 

Let us derive equations of motion for the quantum me- 
chanical amplitude Ak(t) which follow from the action 
(1661). One finds 



A k - jA k + k 2 FA k = 0. 



Introducing the "scale factor" 
1 



b(t) = 



1 if t < 



we can rewrite this equation as 



S^ e t/2R S if t > 



b k 2 
A k + 2-A k + -^A k = 0. 



(71) 



(72) 



(73) 



One can immediately come to an important conclu- 
sion: time- dependent amplitudes of the modes of the field 
f(t,r) in the collapsing CGHS background change as if 
the field was propagating in the fixed dSs background 6 



ds z 



dt l 



2Ht 



(dx 2 



dy 2 ) 



(74) 



with the horizon size 1/H = 2Rs. Although it might 
seem surprising, it is not completely unexpected since 
we know that the "in" modes of the field / are subject 
to the exponential redshift when they pass through the 
collapsing shell. 

As usual, instead of solving the Eq. ([73]) it is conve- 
nient to deal with the equation of motion rewritten in 
terms of the "conformal time" 



r/ 



2\x+R s 

p£/2i? s + l/2' 



(75) 



as t — >• +00. For Uk(r]) = Ak(rj)/ yHrf) one finds in the 
asymptotic future limit 



(76) 



v ]Uk = - 



This equation can be solved exactly in terms of Bessel 
functions. The expression for the field f(t,r) is given at 
arbitrary times t by 

f(t,r) = T] (ckH^ikr,) + d k H[ 2) (k V )) x 

/ /D \ikR s 

x(e r/Rs -e) +c.c, (77) 



where n is given by the Eq. ([75]) . and c&, dk are arbitrary 
constants determined by the initial state of the field. 

It is useful to note that the physical frequency of the 
modes, i.e., the one measured by a future asymptotic 
observer is simply given by 



uo = kn 



^phys • 



(78) 



IV. STANDARD CALCULATION OF 
BOGOLYUBOV COEFFICIENTS 

The manipulations with the action of the scalar field 
/ performed in the previous Section are rather formal. 
In order to check whether they have physical sense, we 
shall now calculate explicitly the Bogolyubov coefficients 
between in- and out-vacuum states and show that the 



6 More precisely, its planar patch. Moreover, since b(t) — > 1 at 
t — > -co, the contracting part of CLS3 is absent. Instead, the ex- 
panding part of dSs is smoothly connected to Minkowski space- 
time at t ~ 0. It is also important to keep in mind that space- 
dependent parts of the modes /fc(r*) for our background space- 
time of course differ from the ones of CIS3 (where fk (r) are plane 
waves) . 
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number of produced particles is given by the Planck dis- 
tribution if we follow the usual considerations p], [l9| . 

At t — »• — oo the spacetime is Minkowski up to expo- 
nentially small contributions to the metric. Therefore, 
we can write for the "in" modes 



u k (t,r) 



1 



-ikt+ikr 



2k 



(79) 



As for the out-vacuum, in the previous Section we have 
convinced ourselves that the time-dependent amplitudes 
Ak (t) change with time as if the scalar field is propagating 
in the dSs background spacetime. Correspondigly, it is 
natural to assume that the out-vacuum of the problem 
coincides with the Bunch-Davies (Euclidean) vacuum of 
the scalar quantum field theory in dSs. Note that this 
concerns only the time-dependent part of the field, while 
the r-dependent part is determined by the modes ([38]) . 
Thus, the "out" modes are 



v k .(t,r) = -y/wHl 2) (kri)(e r ' Rs 



e) 



ik'Rs 



(80) 



where r\ is given by the Eq. ([75]) . 

Since we would like to reproduce the standard an- 
swer for the Bogolyubov coefficients and the spectrum 
of Hawking radiation, we will only be interested at this 
moment in modes with krj ^> 1, the case when geometric 
optics arguments can be applied. One approximately has 
for such modes 



1 



'2k 



/Rs 



ik Rs 



R s^~ t/2Rs 



37ri\ 



• (81) 



Similarly, one has 



Pk'k = —(v k ',u* k ) 



iRs 

47T 



^B(ik'R s + 1, -i(k + k')R s ) x 
x exp \—2ik Rs _i/ 2 e 



L (3 -t/2R s 



i^r — ikt 
4 



. (84) 



It is clear that the relevant 7 parts of the expressions ([83] 
and ([84]) coincide with the ones presented in [l9| . There- 
fore, the analysis of Giddings and Nelson can be com- 
pletely applied to our case revealing 



\Pk'k\ 



e 4Trk'R s — 1 ' 



(85) 



the Planck spectrum of Hawking radiation from CGHS 
collapsing shell. Also, introducing the basis of 
wavepacket modes [19| and tracing out the states local- 
ized in the interior of the shell, one finds that the density 
matrix is completely thermal. 

At this point, we remind ourselves that the two key 
points of this analysis are (a) using the near horizon ap- 
proximation (this also provides a bridge to the canoni- 
cal result [l| for (3+l)-dimensional Schwarzschild black 
holes) and (b) manipulating the contour of integration in 
QLwk and /3k'k i n the complex plane. Namely, noting that 
the main contribution to ak>k and fiyk comes from the 
region near horizon, one finds 



Otk'k 



+oo 



dAre' 



-ikAr ik' R s log (Ar/R s ) 



To calculate the coefficients of the Bogolyubov trans- 
formation between in- and out- vacua, we pick a spacelike 
hypersurface t = Const. ^> Rs and calculate the scalar 
products between "in" and "out" modes on this hyper- 
surface. One finds 



and 



Ztt 



^C dre 



-ikx 



/ ,„ \ ik'Rs 

(e r/Rs - ej x 



x exp ^ik'Rs^e-^ +^+ iJdj . (82) 

The integral in the last line can be taken exactly giving 
iRs fie 



= - —y —B(ik'R s + l,i(k - k')R s )x 
exp f^ik'Rs^e-^ + ^ + ikt) , (83) 



where B(x, y) = ^2+^ ^ s beta-function. 



Pk'k 



r+oo 

I ( iAre ikAr e ikfRslog ( Ar / Rs \ (87) 



where Ar = r — Rs. The arguments of the exponent in 
the integrands of ([86]) and ([87]) seem to be analytic ev- 
erywhere in the complex plane of Ar except the negative 
real axis. Then, we deform the contour of integration in 
([87]) to the negative Ar axis and change the variable of 
integration as Ar — — Ar. This gives us the relation 

\a k , k \^e 2 * k ' Rs \(3 k , k \, (88) 

which together with the normalization condition for the 
Bogolyubov coefficients immediately leads to the Planck 
spectrum ([85]) . 



7 That is, contributing to the absolute values |ajfc/fe| and \fik'k\ °f 
Bogolyubov coefficients. 
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CGHS COLLAPSE: QUANTUM KINETICS 
OF OUT- VACUUM FORMATION 



In the calculations of the previous Section it has been 
assumed that the out- vacuum of the scalar field / coin- 
cides with the Bunch-Davies vacuum of dSs. We would 
like to check now whether this is indeed so, since the 
Heisenberg-like formalism of quantization ([36]) allows us 
to study the very process of formation of the out- vacuum. 

Our program to address this problem is the follow- 
ing. First, in the Sec. IV Al we derive the functional 
Schrodinger equation for the evolution of the quantum 
state ^>(A k ,t) of the scalar field f(t,x). This equation 
describes the first-quantized version of the scalar field 
theory on the background created by the collapsing shell. 
Due to the time-dependence built into the problem, the 
corresponding evolution operator is not diagonal in the 
basis of "in" states. We diagonalize it in the Sec. IV Bl 
and then derive the quantum kinetic equation for the oc- 
cupation number of states instantaneously diagonalizing 
the evolution operator. We shall further denote elemen- 
tary excitations described by these states as "instanta- 
neous excitations". In the Sec. IV CI we exactly solve this 
quantum kinetic equation and show that the out- vacuum 
of the field / is formed as Bose condensate of instanta- 
neous excitations. We then study the properties of the 
out- vacuum. 



A. The Hamiltonian and functional Schrodinger 
equation 

Let us start by deriving functional Schrodinger equa- 
tion for the evolution of the quantum state \£ of the scalar 
field f(t,r) H3. 

Redefining the time variable according to 



the functional Schrodinger equation: 



f dt f F(t f ) = f 
Jo Jo 



dt f 



(89) 



we find from the action ([66]) that the Hamiltonian of the 
field / is given by 



H 



\ J dk(\A k \ f 



where 



F(t) 



2 1 ^l\A k \ 2 ) 



(90) 



(91) 



and prime denotes derivative with respect to the time r. 

The first quantization of the field theory / on the col- 
lapsing CGHS background can be performed by intro- 
ducing the Fock space of states diagonalizing the Hamil- 
tonian (|9U|) . Such states evolve with time according to 



Hip 



1 



d 2 



2d(ReA k ) 2 2 

.diP(A k ,r) 



(ReA k ) 2 



(Re -> Im) 



dr 



(92) 



where A k is as usual the amplitude of the eigenmode with 
momentum k ^> 1/Rs- Since ([90]) is simply the Hamilto- 
nian of a harmonic oscillator with a time-dependent fre- 
quency, the time-dependent quantum mechanical ampli- 
tude A k (r) has a physical meaning of coordinate variable 
for this oscillator, while A k (r) has a physical meaning of 
momentum canonically conjugate to the coordinate vari- 
able. 

Note that because of the time-dependence of the back- 
ground one could introduce an infinite family of different 
Hamiltonians related to each other by time-dependent 
canonical transformations [28|. Among those, we hand- 
pick the Hamiltonian ([90]) because of its simplicity and 
the direct relation to the action ([66]) . As will become 
clear from further discussion, different Hamiltonians re- 
lated by the time-dependent Bogolyubov transformations 
of course correspond to different numbers of elementary 
excitations N k (t) and different histories of N k (t), but 
they provide the very same answer for the final state, 
out-vacuum, of the scalar field /(£, r) at t —> oo. There- 
fore, since our only goal is to determine the structure 
of the out-vacuum state, the choice of H from the class 
of Hamiltonians related to each other by time-dependent 
canonical transformations is largely arbitrary. 



B. Deriving the quantum kinetic equation 

While it is possible to solve this time- dependent 
Schrodinger equation directly [13, , we shall develop 
below a different, more physically transparent approach 
allowing us to find an explicit expression for the occu- 
pation numbers N k . Namely, let us consider Ehrenfest 
equations of motion for the time-dependent operators 
A k (r) following from the Hamiltonian ([90]) : 

—j— = p k , and — = -fi (T)A k . (93) 

Instead of A k and p k let us introduce another pair of 
operators 



1 



(p k - iVt(T)A k ) and 



(p fe + m(T)A k ) 



(94) 



It is clear that the operators b and b^ are in the same 
relation to the operators A k and p k = —id/dA k as the 
creation and annihilation operators to the operators of 
momentum and coordinate in the case of a single quan- 
tum mechanical harmonic oscillator. These b operators 
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diagonalize the instantaneous Hamiltonian ([90]) for all 
values of the time dependent angular frequency f2(r), and 
the combination btf gives the corresponding occupation 
number for a mode with momentum k. 

By explicitly taking derivatives of (|M]h one finds equa- 
tions of motion for the operators b: 



db_ 



- ft' * + 



dr 2fl 



(95) 



Note that we have additional anomalous terms in the 
R.H.S.'s of the Eqs. ([95]) due to the time dependence of 
the frequency ft. 

Now, if we define the occupation number as Nk = bkb\ 

and an anomalous operator xt = Hfi-k (^ ne l^ter de- 
scribes pair creation and is not zero due to the time de- 
pendence of ft), and using ([95]) we come to the three 
coupled first order differential equations 

-^(xl+X fe ) 06) 



dNk 
dr 



^ = -2^-^ + 1) 



^(2^ + 1). 



(97) 



(98) 



where primes denote derivatives with respect to r. In- 
stead of Xk an d Xk ^ * s more convenient to use the op- 
erators Xk and Yk defined according to 

Xk = Xk + xh 

Y k = i( X k ~ x{) (99) 
In terms of them, we can rewrite ([96]) - ([98]) as 



dNk 

dr 2JT2 h 



dXk 
dr 

dr 



2SlY k 



= — l^lXu 



(100) 
(101) 
(102) 



Here we note that ([92]) is written in terms of r, not 
asymptotic observer time t. Hence, the frequency ft is 
not the physical frequency. From ([89]) we see that in 
respect to time of observer living at spatial infinity, we 
can write the physical frequency as 

u>(t) = ky/F® = ± (103) 



Using ([T3T]) and ([T03]) we can also rewrite ([TOO]) - ([T02]) 

in terms of the observer time t as 



Nk = —Xk, 

X k = 2ujY k + - (2N k + 1) 

UJ 

Y k = -2toX k . 



(104) 

(105) 
(106) 



Let us briefly discuss the physics behind the system of 
equations (|104|) - ((106)) . The expectation value 8 of the 
operator Nk is related to the occupation number for the 
mode with a given momentum k. If the frequency uo does 
not depend on time, the system ()104j) - (jlQ6j) reduces to 



N k =0, 
Xfe = 2cuYk, 
Yk = -2ujX k , 



(107) 



i.e., the occupation number remains the same while the 
expectation values of the anomalous operators Xk and 
Yk oscillate with time. This is exactly what one should 
expect, since the time dependence of the Heisenberg op- 
erator bk(t) in the static situation is simply determined 
by the phase factors ex.p(—iHt), and the anomalous op- 
erators Xk and Yk are trivially related to the Heisenberg 
operators bk and b\ (see ([99]) ). When the frequency lj 
becomes time-dependent, spontaneous particle produc- 
tion from the vacuum is possible, with the amplitude of 
the corresponding process being proportional to the adi- 
abaticity parameter uo/uo. This spontaneous particle pro- 
duction is the reason why the anomalous terms appear 
on the right hand side of (fKMl) - ([T06]) . 

The expectation values written in terms of the "in" 
momenta k are of little interest for an observer located 
at spatial infinity who sees the rays passing through the 
collapsing shell, as in the original picture considered by 
Hawking. The frequencies of these rays are exponentially 
strongly redshifted according to (|103|) , so in order to con- 
struct the map between our Heisenberg-like picture and 
the standard calculation of Bogolyubov coefficients, we 
have to rewrite our observables as functions of time t 
measured by a clock at spatial infinity and the physical 
frequency u(t) given bv !103p . 

It is easy to check that the operators Nk , Xk and Yk do 
not commute with each other [30]. What this fact really 
tells us is that for a given observer it is never possible to 
measure simultaneously the occupation number for the 
given mode and the corresponding quantum phase of the 
state. Also, it is impossible to prepare a state with a 
fixed occupation number and a fixed value of the quan- 
tum phase: we can either prepare a Fock state with the 
fixed TV but undetermined phase or a correlated state 
with a fixed value of the quantum phase which is a super- 
position of Fock states with fixed NkS. However, what 
the equations ([104)) - (|106j) want to tell us is that the pairs 
of quanta of "instantaneous excitations" are produced in 
a correlated state (namely, the quantum phases of both 
particles in the pair are correlated), and this correlation 
remains preserved with time. In fact, it is possible to 
exclude Xk and Yk from the system (|104j) - ([106)) writing 



Recall that we use the Heisenberg-like picture where states are 
time-independent, while the operators are time-dependent. 
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a single integro-differential equation for the operator 
It has the form 

^r = W)Jo dt Mf) {2Nk(t) + 1)x 

x cos(2[9(f) - 8(t')]), (108) 



where 



S(t) -e(t') = fdt"oj(t"). 



(109) 



This integro-differential equation is nothing but a 
quantum kinetic equation describing production of in- 
stantaneous excitations on the collapsing background. 
Keeping in mind that it is equivalent to the system of 
three linear differential equations (|104p - (|106p , we are 
now ready to study the solutions of (j!08|) . 



C. "Out" vacuum as a condensate of instantaneous 
excitations 



initial condition for the occupation numbers Nk at t = 0. 
Substituting k = zexp (t/2Rs)/2Rs, one finds the actual 
dependence of Nk on time. 

The vacuum-like initial conditions Nk(t = 0) = for 
the quantum kinetic equation (j!08p can be set as fol- 
lows. In principle, we are only interested in the late-time 
behavior of Nk(t,z) at physical frequencies uo ^> 1/Rs, 
i.e., z ^> 1. Apparently, large physical, "out", momenta 
correspond to even larger "in" momenta since 



K z = 2ujR s = 2kR s e~ t/2Rs <C 2kR s . 



(113) 



Therefore, to fix the values of a& and bk by the initial 
conditions at t —> 0, one can safely use the expansion 
of Hankel functions in (|110p at large values of argument. 
Neglecting rapidly oscillating factors, one then finds at 
t 



N k (t = 0) = AV 



4k 2 R s 



(kl 2 + l^l 2 )-^ (H4) 



Since the Hamiltonian equations of motion for the 
mode amplitude Ak(rj) following from the Hamiltonian 
([90]) can be solved exactly, so can the quantum kinetic 
equation (jlQ8j) for the occupation number Nk(t) of the 
instantaneous modes. If we introduce the new variable 
z = 2(j(t)Rs, the exact solution reads 



N(t,z) 



k 2 R s L 



dA k 



\A k 



1 

2' 



(110) 



where L is the linear size ( "volume" ) of the system 9 which 
we shall omit further, 



A k =a k zH[ 1 \z) + b k zH[ 2 \z) 
(see also Eq. ([77ft ) and 

A k {z) 



a k H[ 1) (z) + b k H^>(z) 



(2), 



(111) 



(112) 



+ I (a k (H^(z) - H?\z)) + b k {H<?\z) - H?\z)) 

The time dependence of the occupation number here 
is implicit: while the constants a& and bk are z- 
independent, they still depend on the value of "in" mo- 
mentum fc, because their values are determined by the 



9 One can understand its origin as follows. The Hamiltonian of 
the field is given by (|90|) . Calculating the energy expectation 
value, one has 



(H) = 5(0) J 



dkQ k N k 



where the delta function appears due to its presence in the rela- 
tion (b^b^,) = N^5(k — k r ). Its dimension is inverse energy and 
the physical meaning is the linear size of the system, as becomes 
manifest when the discrete normalization for time-independent 
modes fkir) is used. 



Ctk 



\h\ 



4k 2 R s 



(115) 



in particular, for the vacuum-like initial condition one 
has 



M 2 + M 2 



8fc 2 R s 



(116) 



As we see, to the leading order in 1/k, the complex co- 
efficients dk and bk cannot be separately fixed by the 
initial condition Nk(t = 0) = 7Vo,/c, only the combination 
Kl 2 + \bk\ 2 can. 

Let us first consider the case of vacuum-like initial 
conditions 7Vo,/c = 0. Apparently, the pair a& = 0, 
bk = \/tt /8k 2 Rs which belongs to the set defined by the 
Eq. (j!16p corresponds to the Bunch-Davies out- vacuum 
of the field f(t,r). Indeed, using ([Hip one finds 



A k =b k H[ 2 \z) = 



^ rr(2) 



2^Hn) 



H^(krj) : 



(117) 



and the "out" modes of the field f(t,r) coincide with 
(|8Q|) . Generally, for arbitrary and bk satisfying (|116|L 
the time-independent asymptotics of the amplitude Ak is 
given by 



A k 



wFj (H[ 2 \kri) + e^H^ikrj) 



y/b(rj)(l - e ak+a l) 



(118) 



where Reak < and b(rj) = —AR 2 s /r] 2 is the "scale 
factor". When the constant parameters ak = a are 

describes Mottola-Allen 
|32| . The general case 



/c-independent, the Eq. 
^(2)-invariant vacua of dSs [3l|, 
of a/e's depending on k corresponds to arbitrary vacua 
of three-dimensional de Sitter space, not necessarily in- 
variant with respect to any subgroups of the symmetry 
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group 50(1, 3) of dSs. Apparently, asymptotic states de- 
veloped from no n- vacuum like initial conditions Nq^ ^ 
belong to Fock spaces built on top of generalized Mottola- 
Allen vacua. 

As follows from (jllOp , the spectrum of instantaneous 
excitations constituting the out-vacuum approaches the 
time-independent asymptotics at t — >> +oo. It behaves as 
a power law at z ^> 1: 

N(w) ~ {2 J Rs)2 + 0((2uR s )- 3 )- (H9) 

The behavior at low frequencies z — »• (in particular, 
the behavior of the zero mode k = 0) seems to be sin- 
gular, but in reality this singularity is cured by the time 
dependence. We shall now show this explicitly. 

qjR s = 0.1, 0.05, 0.025, 0.01 

N 

6: 
5: 
4: 
3: 
2: 
1 : 



2 

FIG. 3: Occupation numbers N of instantaneous excita- 
tions versus t/Rs for various values of the physical frequency 
uj^Rs in the CGHS collapsing setup. Here the curves are 
higher for lower values of uo^R s in accordance with f|103|) . 
After a finite time of the order Rs, the occupation number 
reaches a plateau and remains constant, up to diminishing 
fluctuations. This can be understood as the occupation in- 
creasing exponentially after some finite time, while the phys- 
ical frequency U C> is exponentially redshifted, causing the 
number of particles in a particular mode to become constant 
over time. This corresponds to the formation of out-vacuum 
of the scalar field f(t,r). 

Instead of using exact solution ([110)) it is useful to de- 
termine the approximate behavior of Nk at low frequen- 
cies directly from the quantum kinetic equation. The 
latter is a power tool for analyzing particle production 
and vacuum polarization in QFTs on curved background 
spacetimes, and this will gain one some intuition regard- 
ing behavior of the solutions of ([108)) . Let us start with 
the zero mode uo = 0. Taking the derivative of ()105p and 
noting that 

o)^/o;\ 2 ^ 1 



t/R s = 10, 12 




FIG. 4: The spectrum of instantaneous excitations N versus 
u (t) R s for t/Rs = 10 (blue) and t/Rs = 12 (green) in the 
CGHS collapsing setup. 



(120) 



N ~ I e - t/2R. ^t/2R s _ ^ 2 } (121) 

or as t oo 

X^-h^ 2Rs . (123) 

We conclude that the occupation number of the zero 
mode of instantaneous excitations grows exponentially 
with time. In fact, as clear from the plot, the spec- 
trum Nk remains time dependent for a finite interval of 
frequencies, but this interval shrinks to zero with time. 
Let us now show how exactly it happens and what is 
the uo dependence of the infrared time dependent part of 
the spectrum. Note that although the fixed momentum 
states with k < R^ 1 we calculate the occupation num- 
bers for are not the true eigenstates of the Hamiltonian, 
they do provide the states that a fiducial observer will 
measure at asymptotic infinity. 

The spectrum does appear to become more singular as 
t/Rs increases. This behavior can again be seen from 




at late times t ^> l/Rs, we can write X as 

X 



X 



where we used (|104p . This has the solution 



x ~^ e -t/2 Rs d_ e t/R s y 



10 12 14 



which gives for N 
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(fTQ4|) - ([T05]) . and at u one finds 



27V ^2i? s acje- £ / 2jRs 



2i? s u; 



~ 2 e t/2fl, 

2# s u; 

X =2i? s a^e- £ / 2 ^ - Jl-JI™* 
2R s uj 

« - 7 c */2fl a 



2i? s cj 



(124) 



(125) 



where a and 7 are constants of integration. The depen- 
dence on uo~ l is explicit. 

The expression (|124p is only valid for the interval of 
uj < Rg 1 exp (— t/2Rs). It is very easy to understand 
the behaviour of the IR cutoff cjir separating IR time- 
dependent part of the spectrum from its ultra-violet (UV) 
time- independent part. The uo = is the fastest grow- 
ing mode as all modes with u > have lower occu- 
pation numbers. Comparing behaviour of the station- 
ary distribution N(uj) at low frequencies to the time- 
dependent behaviour of the zero mode, we conclude that 



^IR 



■exp(-t/2R s ). 



VI. QUANTUM KINETICS OF OUT- VACUUM 
FORMATION FOR 3 + 1-DIMENSIONAL 
SCHWARZSCHILD COLLAPSE 

Let us now consider the case of (3 + l)-dimensional 
spherically symmetric collapsing shell. Many of our con- 
clusions applicable to the CGHS collapse hold for this 
case as well. 

To determine the spectrum of instantaneous excita- 
tions and describe the quantum kinetics of out- vacuum 
formation, we first briefly review the set-up and Hamil- 
tonian for the radiation coupled to the background of a 
collapsing massive shell. 

Once again, the action of the scalar field takes on dif- 
ferent forms, one for the interior and one for the exterior 
of the shell. Near the horizon, where most of the particle 
production will occur, the dominant parts of the interior 
and exterior actions are used to determine the overall 
action and the Hamiltonian. 

Following [27] in choosing the real basis of modes f k {r) 
for the massless 10 scalar field / instead of the complex 
basis analogous to ([38]) , ([39]) and neglecting contributions 
to the action exponentially subdominant at t — »• +00, 
one finds similar to the case of CGHS collapse: 

[ dtJ2[~l;A k ,M kk ,A k 

J k,k' 



A k 'N kk fA k , 



(126) 



1 The case of massive scalar field is considered in [29 



where 



drr 2 f k (r)f k > (r), 



N kk > 



hk> = — f 



(127) 



(128) 



Again, matrices M and N are almost diagonal when both 
k,k' ^> Rg 1 , while infrared modes (with kR s < 1) are 
getting strongly mixed since they do not diagonalize the 
Hamiltonian of the field /. The Hamiltonian for a given 
eigenmode B k is 



H 



1 d 2 



-mn 2 (r)B 2 



(129) 



where m is the corresponding eigenvalue of the matrix M, 
and we have defined the "conformal" angular frequency 
as 



w B mB' 



(130) 



and K is the eigenvalue of the matrix N. The conformal 
time is given by 



J dt'B(t' 



where 



and 



B = 1 



Rs 
R(t) 



R(t) ~R S + (Ro - Rs)e 



-t/Rs 



(131) 



(132) 



(133) 



is the position of the shell. This is again just the Hamil- 
tonian for uncoupled harmonic oscillators with time- 
dependent frequencies. 

The spectrum of instantaneous radiation constituting 
the out- vacuum of the scalar field can be found by solving 
the quantum kinetic equation. Performing the analysis 
analogous to the one presented in the previous sections, 
we find that the quantum kinetic equation describing pro- 
duction of quanta of the instantaneous excitations is re- 
duced to the set of three coupled first order differential 
equations ([TM]) - ([T06]) with 



^0 

a{t)'' 



(134) 



where a(t) ~ exp (t/2Rs) is the "scale factor" of dSs 
space, which the modes A k effectively propagate in. Cor- 
respondingly, solution of these equations shows that all 
the conclusions of the previous sections remain the same: 

• the spectrum of instantaneous excitations quickly 
(within the time scale of the order t ~ R^ 1 ) ap- 
proaches the time-independent asymptotics (|1 10p : 
the form of this asymptotics determines the out- 
vacuum state of the scalar field /(£, r), 
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• the form of the most general out- vacuum state is 
given by the Eq. (fTTH]). 

• the UV part ujRs ^> 1 of the spectrum of instanta- 
neous excitations behaves as a power law (j!19|h 

• the deep infrared part ujRs <C exp (—t/2Rs) of the 
spectrum of instantaneous excitations is not sin- 
gular, instead the occupation numbers there keep 
growing exponentially with time; also, at very small 
frequencies the fixed momentum states are not the 
ones diagonalizing the Hamiltonian ([129)) . 



VII. REVIEW OF RESULTS AND DISCUSSION 

In the present paper we study the gravitational col- 
lapse of an infinitely thin shell in (l + l)-dimensional dila- 
ton gravity (see the classic paper |5[) as well as a spher- 
ically symmetric collapse in (3 + l)-dimensional gravity. 
For both setups, the Schwarzschild-like coordinate sys- 
tem is used where full collapse takes infinite time, but 
the shell approaches its apparent horizon exponentially 
quickly. 

Following Hawking [1] (and mostly using the same ap- 
proximations — the limit of geometrical optics, near- 
horizon approximation, analyticity of observables in the 
complex plane of "out" momentum everywhere except 
the negative real axis, etc.), we consider the quantum 
field theory of a massless scalar field on a fixed back- 
ground spacetime corresponding to a collapsing shell. 
Quantum radiation is emitted during such collapse prior 
to the formation of the event horizon. This leads to the 
"polarization" of the in- vacuum of the scalar field /(£, r) 
(i.e. to the creation of a coherent state of excitations of 
the in- vacuum) and to Hawking evaporation of the col- 
lapsing shell. In this paper, we are especially interested in 
understanding the dynamics of the vacuum polarization, 
although the computation of the spectrum of Hawking ra- 
diation is also performed. Our ultimate goal is to better 
understand the nature of the information loss paradox. 
The method which we develop and use to study dynam- 
ics of the vacuum polarization, deriving and solving the 
quantum kinetic equation for particle production during 
collapse, has to our knowledge never been used in this 
context in the literature. 

The main result of our paper is the following: it is 
shown that within a time scale of the order Rs the out- 
vacuum of the field f(t,r) is formed as a Bose-Einstein 
condensate (BEC) of "in" excitations. Such a conden- 
sate is associated with the out-vacuum because it coin- 
cides with the lowest energy state as t — »• + oo. Contrary 
to the standard assumption about the form of the out- 
vacuum, the BEC can become inhomogeneous depending 
on the history of collapse. Furthermore, the informa- 
tion about this collapse history remains encoded in the 
quantum phase (s) of the BEC, which in principle can be 
measurable as we discuss below. 



A. Main points of the paper 

As the way we come to this conclusion is rather tech- 
nical, it is useful to list and discuss the main steps. 

1. Heisenberg-like quantization of the scalar field f(t,r) 
is used instead of the Schrodinger quantiztion (Eq. (f3T|) 
usually found in literature: the field is expanded into 
space-dependent modes /fc(r), while the time-dependent 
amplitudes Ak(t) of the modes are promoted into quan- 
tum mechanical operators (see (|36|) . This allows one to 
explicitly follow the time dependence of observables that 
are the functions of the operators Ak(i), such as the ex- 
pectation value of the stress energy tensor of the field 
/, and its two-point and multiple-point correlation func- 
tions. We do not explicitly separate the Hilbert space H 
of the problem into interior and exterior regions Hi and 
H e - Although the analytic form of the space-dependent 
modes fk is different inside and outside the shell, the 
quantum mechanical amplitude of the mode Ak(t) is a 
single variable. Consequently, during the collapse evolu- 
tion H never becomes the direct product of Hi and H e 
in the slicing that we use. Not even as t — )• +oo. This 
is an additional argument in favor of unitarity preserva- 
tion. Also, note that this Heisenberg quantization picture 
is very much in the spirit of black hole complementar- 
ity conjecture [15]: the time dependent amplitudes Ak(t) 
"know" about both the interior and the exterior of the 
collapsing shell. 

2. The action of the scalar field f(t,r) receives contri- 
butions from both the interior and the exterior of the 
collapsing shell (see ([44]) ). We show that some of the 
contributions to the action can be neglected in the limit 
t — )• +oo compared to others. Roughly speaking, the ki- 
netic contribution to the action of the scalar field as a 
functional of the amplitudes Ak(t) is dominated by the 
region inside the shell, while the gradient contribution 
is dominated by the region outside the shell. Of course, 
when we neglect some pieces of the action compared to 
others, we explicitly violate the locality of the theory. 
Yet, this is not of great concern, as it remains clear how 
to compute local quantities such as the expectation value 
of the stress-energy tensor using the formalism that we 
develop. Also, it was long believed that the ultimate reso- 
lution of the information loss paradox is related to taking 
non-local effects into account (see for example (l3|). We 
strongly feel that the developed formalism captures the 
leading such effect. 

3. Once only leading (in the late-time limit t — »• +oo) 
terms are retained, the action is shown to be equivalent 
to the action of a massless scalar field on a dSs back- 
ground (see (fT3|) and following discussion). This is not 
terribly surprising. 11 In the limit t — »• + oo, the collaps- 



What perhaps is surprising is the effective dimensionality of this 
de Sitter space; note that the original collapse problem is two- 
dimensional. 



17 



ing shell is located exponentially close to its horizon, and 
the modes of the scalar field f(t,r) passing through the 
shell experience exponential redshift. Since the shell is 
in motion, the amount of redshift itself depends on time. 
It should therefore be possible to quantify the effect of 
redshift by introducing the scale factor b(t) relating the 
frequencies of "in" and "out" modes. The developed for- 
malism allows us to do exactly that. 

4. As spacetime is flat at t — >> — oo, up to exponen- 
tially small corrections to the Minkowski metric in the 
case of CGHS collapse, we assume that at t = the 
quantum state of the scalar field remains in the "in" 
(Minkowski) vacuum for modes with momenta k > Ro 
(where Ro = R(t = 0) is the radius of the shell at t = 0). 
In a sense, the time-dependent amplitudes Ak(t) of the 
modes effectively propagate in a spacetime with scale fac- 
tor 6(t), which behaves exponentially at t — >> -f oo and ap- 
proaches 1 as t — >• — oo, so the effective spacetime is not 
dS%, but the expanding planar patch of dSs smoothly 
connected to the flat Minkowski space at t ~ 0. We 
correspondingly set vacuum-like initial conditions for the 
modes with k ^> Rq 1 of quantum state of the scalar field 
at t = 0. 

5. Assuming that the out- vacuum of the field / is effec- 
tively the Bunch-Davies ("Euclidean") vacuum for the 
scalar QFT on dSs background, we compute the Bo- 
golyubov coefficients between "in" and "out" vacua in 
Sec. HVl for the CGHS collapse and show that they pre- 
cisely coincide with the standard result by Giddings and 
Nelson [l9|, and the Hawking radiation is expected to 
have the Planck spectrum. It is also clear from the 
derivation that the same result applies for the case of 
(3 + l)-dimensional spherically symmetric collapse of a 
thin shell. 

6. Instead of assuming the explicit form of the out- 
vacuum, we would like to derive this form explicitly as 
a result of the evolution of the in-vacuum polarization. 
To describe this evolution, we use the formalism of the 
quantum kinetic equation developed in Sec. IV Bl This 
formalism was first introduced in [33| and mostly used to 
study particle production in the cosmological context (see 
for example recent papers (Hf). Historically, the formal- 
ism of the quantum kinetic equation fell out of favor due 
to the difficulties in prescribing the physical content to 
the definition of particles for QFTs in curved spacetimes 
and likewise to the procedure of instantaneous diagonal- 
ization of the Hamiltonian [28]. In the next subsection we 
explain how to deal with these difficulties and understand 
them, and how to extract the physical content from the 
quantum kinetic equation describing particle production 
in a QFT on a fixed curved background. 

7. Solving the quantum kinetic equation exactly in Sec. 
IV CI we find that, within a time scale t ~ Rs (counted 
from the initial moment of time t = 0), the Bose-Einstein 



condensate of "in" excitations 12 is formed. This conden- 
sate can be interpreted as the out- vacuum as it is the low- 
est energy state of the field at t — )• + oo. The form of the 
condensate is given by Eq. ([11 8|) . For general conditions 
at t = 0, it corresponds to generalized Mottola- Allen 
vacua of three-dimensional de Sitter space, the Bunch- 
Davies ("Euclidean") vacuum being just one particular 
case. The word "generalized" used here means that the 
complex quantity a parametrizing Mottola- Allen vacua 
can be an arbitrary function of in- moment a k. 

We shall now discuss the implications of some of the 
observations listed above. 



B. Notes on the quantum kinetic equation 

A quantum kinetic equation like (|108|) is an equation 
driving the evolution of occupation numbers Nk(i) asso- 
ciated with a Hamiltonian H (roughly, H ~ J dk ■ cj/cA^, 
where uok is the physical, observed frequency of excita- 
tions). As the Hamiltonian H can only be defined in 
QFTs on fixed curved backgrounds up to time-dependent 
canonical/Bogolyubov transformations of the field (28| . 
the occupation number by itself is not a physical, 
observable quantity. However, it is directly related to 
observables such as the expectation value of the stress- 
energy tensor of the field / or its two-point correlation 
function. For example, in our setup the (£, t)-component 
of the stress-energy tensor is given at late times by 

(Ttt) = j cLu k oj k [(2N(t,co k ) + l)/ fc (r)/ fc *(r)+ 

+ X(t,u)f k (r)r_ k (r)] + 

+ ^L= J dco k Lo k [(2N(t,co k ) + l)/ fe (r)/ fc *(r)- 
-X(t,u)f k (r)r_ k (r)], (135) 

while the two-point function of the scalar field / at equal 
(late) times is 

{f(t,r)f(t,r')) = 

-X k (t,w k )f k (r)f*_ k (r')], (136) 

where N and X are the solutions of Eqs. ([104)) and (|105|) . 

Although both expressions ()135|) and (j!36p are not reg- 
ularized, it is often useful to keep UV divergent terms in 
them explicitly (for example, in order to see how UV di- 
vergences build up with time). To regularize expressions 
(j!35p and (jl36j) . one either has to introduce a physical 



Or, more precisely, quanta of states instantaneously diagonaliz- 
ing the Hamiltonian of the scalar field /. 
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UV cutoff or introduce the adiabatic basis of eigenfunc- 
tions/ fc Blilil. 

Quantum kinetic equations similar to (|108|) exist (and 
their structure is nearly universal) for many physical ef- 
fects involving time-dependence of a background field. 
(See, for example, [35j, where a similar quantum kinetic 
equation is derived for the process of Schwinger pair pro- 
duction in a strong electromagnetic field.) In particular, 
such an equation can be derived for an arbitrary quan- 
tum field theory in an arbitrary background spacetime. 
If a theory does possess an out- vacuum/ vacua, described 
by a quantum kinetic equation similar to ([108)) , then the 
equation admits a time-independent solution(s). Such 
solutions are attractors in the sense that general solu- 
tions of (|108p approach them as t — >> +00. The physical 
meaning of such attractor solutions is that they encode 
the effects of in-vacuum polarization (i.e. the develop- 
ment of a new vacuum, different from the "in" state) 
driven by the time-dependent background field. 

In principle, since Eq. ([108)) completely describes the 
physics of particle production in a background spacetime, 
its solutions also contain information about particle pro- 
duction, the Hawking effect in our particular case. 13 The 
Hawking effect is encoded in time-dependent corrections 
to the time-independent asymptotic solution of (|108|) . 
These corrections to the numbers Nk(t) of instantaneous 
excitations vanish as t — )• +00, but their contribution to 
the expectation value of the energy- momentum tensor re- 
mains finite. For example, as can be seen from ([135p , cor- 
rections to the time- independent spectrum N(uj), which 
scale with time like exp (—t/2Rs), lead to a finite energy 
flux at spatial infinity and describe the Hawking effect. 
This observation also makes it manifest that the Hawk- 
ing effect is subdominant with respect to the effect of 
Bose-Einstein condensation leading to the formation of 
out- vacuum. 



C. Properties of the out- vacuum and observables 

Why can the out-vacuum ([118)) be associated with a 
Bose-Einstein condensate of instantaneous excitations? 
In the context of condensed matter physics, one usually 
says that the Bose-Einstein condensate is created if the 
spectrum of bosonic (quasi-)particles acquires the form 

n(k) = n S(k) + — , (137) 

CT — 1 

that is, when a sharp peak is developed at zero frequency 
and the occupation number of the zero mode is extremely 



large. We instead have a smooth distribution which be- 
haves as a power law (j!19p at large frequencies and re- 
mains finite (but large) at uok — 0. 

However, in real physical systems the distribution 
(P~37j) is almost never seen in its canonical form. If one 
starts from an arbitrary non-equilibrium initial distru- 
bution n(t,k) (or just cools a system down), a so-called 
"precondensate" is formed first [36| within the charac- 
teristic (short) time scale of the order of the mean free 
time in the interacting Bose gas describing the quasi- 
particles of the system. The precondensate is charac- 
terized by large (but finite) occupation numbers n(k) of 
quasi-particles in the IR part of the spectrum, but the 
distribution n(k) is non-singular. Instead, the behavior 
of n(fc) is smooth at k <C T. 

Another characteristic feature of the precondensate is 
a strong quantum coherence among the states of quasi- 
particles constituting it. However, the global quantum 
phase of the BEC is not yet formed, and the time of its 
formation is much larger than the mean free time and 
actually depends on the system size [36j. To conclude, 
for a physical system with infinite 3- volume the quantum 
state describing a precondensate is as close to the actual 
BEC as one can get even as t — )• +00. 

We would like now to discuss thermal properties of the 
"out" vacua. The Bogolyubov coefficients between the 
"in" (Minkowski) vacuum and the Bunch-Davies vacuum 
(|117|) were calculated in Sec lIV) and it was found that the 
corresponding spectrum of Hawking radiation is thermal 
with the correct temperature. Let us also calculate the 
detector response function for a fiducial observer located 
at r ^> Rs. For a detector interacting with the field 
according to 

J dtm(t)f(t,r(t)) (138) 

(where r = r(t) is the trajectory of the detector, and for 
us is a constant), one has [H, [39| for the probability of 
transition between the states \Ei) and \Ej) of the detector 

/+00 
dte- iAEt G(t,r(t);0,r(0)), (139) 
-CO 

where G(x, y) is the Wightman function of the field 
f(t,r) and = (Ei\m(0)\Ej) . For the Bunch-Davies 
vacuum, the Wightman function G is given by 14 

G= ( +CC dk^H[ 2 >(k V )Hi 2 \k V '), (140) 



As usual, the fundamental difference between vacuum polariza- 
tion and particle production effects is that the former vanishes 
when the background field is turned off (the out-vacuum coin- 
cides with the usual Minkowski 0(d — 1, l)-invariant vacuum), 
while the latter remains. 



It is worth emphasizing again that this G differs from the Wight- 
man function of the scalar field f(t, r) in dS3 : because, for a given 
mode with momentum k, the phase volume spanned by the mode 
is different in our collapse problem and in (IS3. This is why we 
must compute the Wightman function explicitly instead of sim- 
ply borrowing the known result for CIS3. 
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where rj = —2Rs, v' — —2Rsexp(—t/2Rs) are conformal 
times. Note that the lower limit of integration is taken 
to be R^ 1 , as the modes with lower momenta are out 
of equilibrium and are not in the Bunch-Davies vacuum 
state. The resulting expression (j!40p is periodic in imag- 
inary time t with the period AttRs- Correspondingly, the 
detailed balance equation for the detector [38] is 



P 



exp (-ATrRs\Ei-Ej\) = exp (-4irR s AE), (141) 



and we conclude that the detector response is thermal 
with T = Tjj 



AttRs ' 

This is no longer so when the out- vacuum of the field 
/(£, r) is the generalized Mottola- Allen vacuum (|118p . 
Even if a& = a = Const., the detector response is not 
thermal [37], [38|. Indeed, the Wightman function of the 
field / is now given by 

Ga(x,x') = - _ 1_ (G E (x,x') + e a+a *G E (x',x)+ 

+e a 'G E (x, x' A ) + e a G E (x A ,x')), 

(142) 

where Ge is the Bunch-Davies Wightman function (114Q|) 
and the index A means that one takes a point on dSs 
antipodal to x. Correspondingly, the detailed balance 
equation is 



P 

r ''■'.) _ e -4irR s AE 
P 



1 



a+2irR s AE 



1 



ct-2irR s AE 



(143) 



i.e., the detector does not equilibrate, and the Hawking 
flux associated with the vacuum (|118|) is not thermal. 
This observation has important implications for the in- 
formation loss paradox as we discuss in Section IVII Dl 
below. 

Given ([143)) , the probability to excite states with Ej ^> 
Rg 1 is finite for generalized Mottola- Allen vacua |37j . 
which means that the UV tail of the distribution of Hawk- 
ing radiation is much heavier for such vacua than for the 
Bunch-Davies vacuum. 15 



15 It is possible to see immediately, without explicitly solving the 
quantum kinetic equation, that the spectrum N(qjj~) should be- 
have as a power law at large frequencies. The equation for the 
amplitude Ak(t) of the mode, rewritten in terms of the variable 
r defined in (|89J) , has the form of the equation of motion for a 
harmonic oscillator with the variable frequency fifc(r) denned in 
(EJ: 

A'£ + n 2 k (T)A k = 0. (144) 

The occupation number (r) of instantaneous excitations coin- 
cides with the adiabatic invariant associated with this oscillator. 
For large physical frequencies uo^ ^> Kg 1 ■> when the adiabatic 
approximation holds very well, the change in the value of the 
adiabatic invariant AN k (roughly, the number of particles pro- 
duced) is determined by the singularities of fi|(r) (and ^'(t)) 
in the complex plane of r [4(|. If ^K r ) does not have singular- 



D. Implications for the information loss paradox 

The emerging picture sheds some new light on the 
information loss paradox and proposed ways to re- 
solve it. We claim that, for both two-dimensional 
CGHS collapse and (3 + l)-dimensional Schwarzschild 
spherically-symmetric collapse, the time-dependent am- 
plitudes Ak(t) of the modes of the scalar field f(t,r) 
change with time as if the modes were effectively propa- 
gating in a dSs spacetime. The Hubble parameter of this 
effective spacetime is related to the horizon scale in the 
collapsing geometry 



H 



1 



2R S 



(145) 



In the usual quasi-de Sitter geometry describing an inflat- 
ing universe, the modes are constantly leaving the hori- 
zon. At the same time, the value of H keeps falling until 
the slow roll condition 



H 2 



« 1 



(146) 



ceases to hold, and inflation stops. During the subse- 
quent era of power law Hubble expansion inflationary 
modes reenter the horizon and become the seeds for pri- 
mordial perturbations of the matter density Sp/p. If the 
power law Hubble expansion continues ad infinitum, all 
inflationary modes eventually reenter the horizon, so that 
the information about the initial state for inflation be- 
comes completely recovered at t — )• +00. The information 
loss paradox therefore does not appear in this situation. 

What happens in the effective dSs space associated 
with the dynamics of modes on the collapsing background 
once the backreaction effects are taken into account? As 
the collapsing shell/black hole Hawking evaporates, the 
horizon scale Rs decreases. But then the associated 
H scale (|145|) increases] This setup looks like inflation 
driven by an inflaton field <p that rolls up its potential 
V((j)) ~ H 2 = (2Rs)~ 2 towards the super-Planckian 
scales instead of rolling down towards the minimum of 
V ((/)). Correspondingly, the modes will keep leaving the 
horizon even more effectively at late times and will never 
reenter it again. 

Despite this difficulty, we feel that our picture might 
provide an effective resolution of the information loss 



ities on the real axis, is exponentially small since its value 

is given by the singularity closest to the real axis. However, in 
our case ^|(t) ~ k 2 /r, so that there is a simple pole present at 
r = 0. This pole corresponds to the location of the event horizon. 
In this case, ANf~ decays only as a power law at large k. It is in- 
teresting to note that if the horizon is not ultimately formed (for 
example, due to backreaction effects), then the out-vacuum is 
Minkowski, the angular frequency f2fc(r) has no singularities on 
the real r axis, and we return to the exponentially small answer 
for ANk, reproducing Hawking's answer already at this simple 
level. 
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paradox. To see how, let us discuss again the structure 
of the generalized Mottola- Allen vacua ([118)) . The late 
time general asymptotic solution of the quantum kinetic 
equation ([108)) is given by ([118)) . Setting vacuum-like 
initial conditions at t = 0, one can fix only the combi- 
nation \ak\ 2 + |frfc| 2 , not the complex coefficients a& and 
bk- In order to determine the values of a& and bk sep- 
arately, that is to find which particular vacuum among 
the generalized Mottola- Allen vacua ([118)) is realized at 
t — »• +oo, one must fix the 2-point correlation function 
of the field f(t, r) together with N k (t = 0) at t = 0. Such 
correlation functions do not have to respect Minkowski, 
E(2) or any other symmetries at t — »• as they are only 
determined by the initial conditions for the field. 

Imagine now that we initially have a thin collapsing 
shell with radius R(t = 0) = Ro. The shell is allowed to 
have wiggles/internal structure, with information about 
them being effectively encoded at t = in the modes of 
the field /(£, r) with k <C Rq 1 . Note that the state of the 
field will not be vaccum-like for such modes. Let us also 
assume that a correlated EPR pair of quantum particles 
is initially placed right behind the shell, say, at R = 
Ro + SR. Information about the pair and its quantum 
correlations is encoded in the two-point function of the 
field /(£, x), namely, in the amplitudes and the quantum 
phases of the modes with k ~ / _1 <C Rq 1 . 

As we argue, once the collapse has started, the con- 
densate of "in" particles (i.e, the out-vacuum) builds up 
during the time interval of the order ~ R$. Note that the 
shell is still very far from its own apparent horizon at this 
moment. The out- vacuum is the generalized Mottola- 
Allen vacuum (|118[) parametrized by the complex quan- 
tity a k trivially related to the complex coefficents a& and 
bk in ([111)) and therefore encoding the information about 
wiggles/internal structure of the collapsing shell as well 
as the coherent pair of quanta following the shell. 

As was explained in the previous Section, such a vac- 
uum emits a non-thermal flux of Hawking quanta, which 
can therefore carry away information about The UV 
tail of the distribution of emitted quanta is heavy (37| . 
so we expect this vacuum to evaporate away in finite 
time. Since it takes an infinite amount of time for the 
collapse to finalize in the slicing that we use, all the in- 
formation about the wiggles/internal structure/quantum 
correlations of the shell as well as quantum correlations 
of the pair following the shell will be released to spatial 
infinity before the event horizon/black hole is formed. 

Finally, the only information which survives at t —> 
+ oo without getting released to spatial infinity will be 
the information about the total mass of the shell (in 
the case of (3 + l)-dimensional Schwarzschild collapse) 
or the value of the cosmological constant (in the case of 
(1 + l)-dimensional CGHS collapse). The associated out- 
vacuum will have the Bunch-Davies form ([117)) leading 
to the thermal flux of Hawking radiation from the shell. 

The AdS / CFT correspondence in a certain sense seems 
to confirm this scenario. According to it, the process of 
gravitational collapse on an AdS background corresponds 



to the process of thermalization in the non-equilibrium 
N = 4 Super- Yang-Mills (SYM) quark-gluon plasma [4~0 - 
13 • The N = 4 SYM theory is a conformal field theory 
(CFT), and CFTs are known to take a very long time to 
t her malize starting from a non-equilibrium state. (See for 
example j45[, where thermalization of a classical scalar 
CFT is studied.) Typically, when one starts from an ar- 
bitrary non-equilibrium state in a CFT and follows the 
kinetics of equilibration, the spectrum of CFT excita- 
tions quickly (within the mean free time for the plasma 
of CFT excitations) approaches time-independent non- 
thermal asymptotics. Such asymptotics are characterized 
by power-law specta /c -7 of bosonic excitations, where 
the exponents 7 are determined by collision integrals of 
the kinetic equation describing such thermalization and, 
correspondingly, interactions of the CFT. Such asymp- 
totic solutions have a very long life time (much longer 
than the mean free time), although the CFT of course ul- 
timately thermalizes. This phenomenon has been called 
"prethermalization" in cosmology and in the heavy-ion- 
collisions literature. 

We suspect that the prethermalization phenomenon is 
exactly what is seen in [4ll-[44{ . The theory quickly (dur- 
ing one mean free time) reaches dynamic equilibration. 
This regime is characterized by components of the stress- 
energy tensor reaching their equilibrium asymptotic val- 
ues, but continued non- thermal behavior of correlation 
functions. 

To conclude, the scenario presented above by itself 
does not probably constitute the resolution of the in- 
formation loss paradox as many steps which we discuss 
should be significantly better quantified. In particular, 
our understanding of the decay of Mottola- Allen vacua 
(the decay rate, the spectrum of decay products, etc.) 
should advance well beyond [37l, the analysis similar to 
the one performed in the present paper should be repro- 
duced for different slicings since the Schwarzschild coor- 
dinate system does not cover the whole BH spacetime 16 . 

Our discussion only involves collapsing shells made of 
the very same matter content as the most rapidly emitted 
modes of Hawking radiation (massless scalars). It would 
be instructive to understand what happens in the case of 
collapsing shells made of dust-like matter (for example, 
massive scalar with m > Rg 1 ) and where exactly the 
information about their internal structure is encoded in 
Hawking radiation. 17 



Especially, it is necessary to determine what an infalling observer 
experiences. 

A preliminary analysis [29| shows that the structure of M h y 
and N k y kernels (see Section [TTT]) for massive scalars is not very 
different from the one for massless scalars (the case discussed 
in this paper). Therefore, we expect that the quantum kinetic 
equation for the number of instantaneous excitations will have 
the form analogous to (|108|) . and our conclusions concerning the 
structure of the out-vacuum will be also similar. Then, the in- 
formation about the initial state of the shell made of massive 
scalar particles will be encoded in the parameters of the gener- 
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Ultimately, one must understand the behavior and ori- 
gin of black-hole entropy using the formalism developed 
here. Nevertheless, we expect that this work is an im- 
portant step towards the ultimate resolution of the cele- 
brated information loss paradox. 
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Appendix A: Massive shell collapsing on a 
primordial black hole 

In this Appendix we first find the classical equations of 
motion for a massive shell collapsing onto a pre-existing 
primordial black hole (PBH). After which, we set up the 
equations for a scalar field coupled to the background of 
the collapsing shell around the PBH. 



1. Classical solution 

Here we take that the internal metric is that of a PHB, 
so that the metric is just described by the Schwarzschild 
metric 



ds 2 



1 



2GM 1 



dT 2 



dr 2 



1 - 2GMi/r 



r W 



(Al) 

where Mi is the mass of the PHB. The exterior metric is 
also given by the Schwarzschild metric, however the mass 
here is now the total mass, i.e. that of the PHB and the 
shell (with mass M2): 



ds\ = 



1 



2GM t 



tot 



dT 2 



dr 2 



1 - 2GM tot /r 



(A2) 

where M tot = M x + M 2 - 

For a given interior metric with metric coefficient 
J(r = R) and exterior metric with metric coefficient 
P(r = R), one can find the perpendicular components 
of the acceleration on either side of the shell. Taking the 
ratio of the sum and difference of these once can find 



4naG(a + /3) 



2(a 2 - p 2 ) 
R 



J' - P' 



(A3) 



alized Mottola- Allen vacuum realized at t — > +00. The associ- 
ated Hawking radiation will be strongly non-thermal, and the 
information about the initial state of the system will be released 
before the collapse is finalized. 



where a is the energy density of the domain wall and 



a = y/J + R 2 and (3 = V J + R 2 . (A4) 
From (lAlj) and (fA2|) we have that 
2GMi 



J = 1 



R ' 



(A5) 



2GM tot _ 1 2G(M 1 +M 2 ) ^ (Ag) 



R 



R 



Substituting (|A5[) and (|A6|) into (|A3[) and solving for the 
mass of the shell gives: 



= ^^(^jTrJ+^/p + rJ 



or solving explicitly for the mass 



(A7) 



M 2 = fsfiZ™f+x> - 2x.GR) . (A8) 

A mass of this form is actually a constant of motion, 
i.e. dM2/dr = 0. Therefore, we can take M2 = H to be 
the Hamiltonian of the system. Note that the Hamilto- 
nian is in terms of the infalling observer time r, an ob- 
server who is riding on the shell as it is collapsing. What 
we want, for Cosmological purposes, is the Hamiltonian 
in terms of the asymptotic observer time t. To get this, 
we must first find the Lagrangian of the system to deter- 
mine the action, then perform a coordinate transforma- 
tion on the action, since the Hamiltonian is not invariant 
under coordinate transformations, and then find the new 
Hamiltonian. To find the coordinate transformation we 
use (|A2j) and find that the relationship between r and t 
is given by 



(A9) 



The Hamiltonian as observed by the asymptotic observer 
is then given by: 



H = WP3/2 J JP2 - {J - P)A2 - 2,oGR 
9VP 2 - R 2 \ V P 2 -R 2 

(A10) 

We are mostly interested in the near horizon limit, R —> 
Rs = 2GM tot , so that in this limit J — >• const and P — >• 0. 
Taking the near horizon limit of ([A 10)) we then have 



H = 



(All) 



9Jf 2 -R 2 



where 



a = a (%y/j(R = R s ) - 2iraGR s ^ 

( I 2GM 1 _ N 

a I 3a/ 1 - 2ircrGRs 



(A12) 
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To find the radius of the shell as a function of asymp- 
totic observer time £, we need to solve (|Alip for which 
gives 



R 4 



(A13) 



where h = 9H/ (4ir/i). In the near horizon limit we can 
then write 



(A14) 



Therefore the dynamics of shell in the near horizon 
regime can be obtained by solving R = — P, where we 
took the minus sign since the shell is collapsing, which is 
identical to that in 1271. 



2. Quantum radiation of instantaneous excitations 

Here we follow the procedure found in [27], where we 
split the action in to two parts, one for the interior and 
the other for the exterior, and find the dominant pieces. 
For the interior metric we have 



2tt / dt 



Rs 



dr r z 



1 



Tl 
f I 



2GM 1 /r 
2GM 1 



(A15) 



and the exterior we have 



S+ = 2ir dt 



dr r 



Rs 



2GM tot /r 
2GM tot 



(A16) 

To move on, we must determine T = (dT / dr)(dr / dt). 
From the structure of the metric, the relationship be- 
tween r and T is analogous to that of r and t (|A9|) , 
which is given by 

^=jVTTB}. (A17) 
(jAlTj) and (|A13|l we find that 
P 



Therefore using (jA 
1 



T = 



J 



JP 



J-P 



R 2 = 



J 



1 + ( J - p )^ (A18) 



so that in the near horizon limit f — >• P '/ J ~ P ', since in 
this limit J is finite. Comparing terms and keeping only 
the dominant contributions, we can then write the total 
action as 

r Rs 



2tt / dt 



H 
L 



dr r 



drr [ 1 



1 - 2GM 1 /r 

- (d r <S>? 
r 



(d r $) 2 (A19) 



As in [27] one can expand $ in modes and define matri- 
ces M and N that are independent of R(t), here defined 

by 

rRs 



M kk , = 4tt 



drr 2 fk{r)fk>{r) 



1 - 2GMi/r' 



(A20) 



Rs 



Rs 



N fefe , = 4tt / dr r 2 1 - -2- f k {r)f k , (r). (A21) 



Using the standard quantization procedure and perform- 
ing a principle axis transformation, we find that for a sin- 
gle eigenmode, the Schrodinger equation takes the form 



J_cP_ 
2m db 2 



m 



Hb,v) = i^ (A22) 
or] 



where fa denotes the eigenvalue of M, b is the eigenmode, 



K 



1 



r,2 



m 1 - R s /R 1 - R s /R 
where K denotes the eigenvalue of N and 



T] 



1-^ 



(A23) 



(A24) 



This has identical structure to that found for the massive 
shell collapsing upon itself. 



Appendix B: Collapse of two shells 

In this Appendix we will again first find the classi- 
cal equations of motion for two massive shells collapsing. 
After which, we set up the equations for a scalar field 
coupled to the background of the collapsing shells. 



1. Classical Solution 

Consider two concentric shells collapsing under grav- 
itation, where the inner shell is located at a radius R\ 
with mass M\ and the outer shell is located at radius R2 
(where R2 > Ri) with mass M2. Thus we have three 
regions denoted by /, //, and see Figure [5] Due to 
Birkhoff 's theorem, the metric in region / is taken to be 
flat Minkowski 

ds] = -dT 2 + dr 2 + r 2 dtt 2 , (Bl) 

the metric in region is taken to be Schwarzschild 



ds 2 n = _ ( 1 _ x*M± \ d p 



dr 2 



1 - 2GM l /r 



r 2 + dti 2 



(B2) 

and the metric in region III is taken to also be 
Schwarzschild 



ds 2 



in 



2GM tt 



dt 1 



dr 2 



1 - 2GM tot /r 



-r 2 +dQ 2 , 
(B3) 
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FIG. 5: Schematic of the collapse of two shells. The interior of 
shell 1, with mass M\ at radius R\ is taken to be Minkowski, 
the region between the two shells is taken to be Schwarzschild 
and the exterior of shell 2, with mass M2 at radius R2 is taken 
to be Schwarzschild. 



where M tot = Mi + M2. Due to the structure of the 
spacetime, each of the shells satisfy the Gauss-Codazzi 
equation separately. Thus considering the shells indi- 
vidually we find that mass Mi, which is a constant of 
motion, is given by 



2tt.5-G.Ri 



(B4) 



where f is the proper time of an observer riding on the 
shell. This result is not unexpected since it is just the 
case considered by Ipser and Sikivie in [26|. Here the 
internal time T is related to the proper time f by 



(B5) 




and the "external" time i is related to f by 




where 



Si EE 1 



2GM X 
Ri ' 



(B6) 



(B7) 



We find that mass M2, which is also a constant of motion, 
is given by 



^ = ^W* 2 +(^) 2 -2naG Rl \ (B8) 



where r is the proper time of an observer riding on the 
shell and 



Bo 



2GM 1 



(B9) 



Again this result is not unexpected since this is the sce- 
nario examined in Appendix [Aj see (|A8[) . The "internal" 
time t and external time t are related to the propertime 
r by 



(BIO) 

(Bll) 

where 




B = 1 



2GM t , 
R2 



(B12) 



Since the masses are separately conserved, we can take 
them as the Hamiltonians of the shells. We are interested 
in the collapse as observed by the asymptotic observer, 
hence we must find the equation of motion of the two 
shells as a function of time t. We will first consider the 
outer shell, M2. 

It is easy to see that there is no difference in the eval- 
uation of the equation of motion as a function of time 
t in the two shell case than that of the shell collapsing 
onto the PBH. Therefore we find that in the near hori- 
zon limit, where here that corresponds to B — >> 0, the 
Hamiltonian corresponds to 



H 2 



4^/ii?^3/2 



Rn 



(B13) 



which is just (jAll|) with the appropriate substitutions 
and fi is defined in (jA12|) . To leading order in £>, the 
equation of motion for the position of the shell is then 
given by 



R2 



-B 



(B14) 



Now we consider the equation of motion for the inner 
shell of mass Mi. The Hamiltonian is a function of f , 
where we want as a function of t. From (|B6|) we see that 
f is a related to t, while from (|B10|) and 



(iBTTj) t is re- 
lated to t via r. Therefore we can perform the coordinate 
transformation 



df df dt dr 
dt dt dr dt 



(B15) 



to transform (|B4|) from a function of f to a function of t. 
By first finding the corresponding Lagrangian and deter- 
mining the action, then performing the coordinate trans- 
formation given in (|B15|h we find that the Hamiltonian 
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as a function of the asymptotic observer time t is given 

by 



B 2 B 2 - (B 2 - B)R 2 2 



B ^ Jb 2 B 2 B 2 - BB\R\ - (B 2 - B)R 



where the superscripts in parenthesis mean the initial 
position of the corresponding mass. Hence we see that 
as t — > oo, i?i — > Rg + const, so that B\ never quite 
reaches zero, however it can be arbitrarily close. 



B 2 B\B 2 — (1 — B^BB^R 2 — B\(B 2 — B)R3^ 2 - Quantum radiation of instantaneous excitations 



B 2 B 2 B 2 - BB\R\ - B 2 (B 2 - B)R\ 
-2iraGR i y (B16) 

We will again be concerned with the near horizon limit, 
where here that corresponds to B\ — ^ 0, which is inside of 
the actual Schwarzschild radius Rs = 2GM tot . However, 
until the Schwarzschild radius is formed, the asymptotic 
observer will be able to see this region. Taking this limit, 
we see that for the squareroot in the parenthesis, the 
numerator and denominator are the same, so (|B16p can 
be written as 



^kgR\B{ 



3/2 



B 2 B 2 



(B 2 



B)Rl 



B ^ \Jb 2 B\B 2 - BB 2 R\ - (B 2 - B)R? 2 
x (1 -2ndGR 1 ) 

Att/IRIB^ 2 B 2 B 2 - (B 2 - B)R 2 



b 2 Vb 



where 



B 2 B 2 B 2 - BB 2 R\ - (B 2 



fl = a(l- 2naGR 1 s ) 



B)R 2 
(B17) 



(B18) 

and R\ = 2GM\. Here we see that (|B17|) depends ex- 
plicitly on R 2 . To simplify things further we will make 
the assumption that R\ and R 2 don't start that far away 
from each other so that we can treat that R 2 is simulta- 
neously in the near horizon limit. Thus we take that R 2 
satisfies (|B14[) that that (|B17[) can be written as 



Solving for R\ yields 
Ri = - 



^liR\B*l 2 B{ 



3/2 



(B19) 



B 2 JB 2 B 2 -B 2 R 2 




BjRfB 3 
B\h\ 



(B20) 



where hi = Hi/^nfi). Since B\ 4 we can Taylor 
expand and see that the dynamics of R\ are given by 

B 



B 2 



Solving (|B2ip to lowest order in Rq 
using the explicit time dependence of B, 



(B21) 
R^ we find, 



Ri ttR^+^'-R^exp 



RsjR^ -Rs) 

H s H s 



( e -t/Rs _ 1} 
(B22) ' 



In this case, the action for the scalar field coupled to 
the collapse of the two shells is then split into three pieces. 
We again must find the dominate parts. The individual 
actions then takes the form 



57 = 2tt / dt 



dr r z 



S n = 2tt / dt J drr 2 



l(d t ®) 2 + T(d r ®) 2 



m) 2 



ll-2GM 1 /r 
2GM 1 



■til 



and 



Sin = 2tt 



dt I drr 2 

R2 



my 



(d r ®f 



(B23) 

(d r ®f 

(B24) 



where we transformed to the asymptotic observer time t. 
Using (fB14|) and (fB14|) we find that f « BB 1 /B 2 and 

t « B/B 2 . Therefore we can write 



2GM tot /r 
2GM tot 



St = 2tt 



J dt drr 2 



S n = 2tt / dt 



R 2 



dr r 



B 2 
BB\ 



Bo 



(d t $) 2 



(dt®) 2 



BB 1 



(dr®f 



(B25) 



B 1 
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B~2 



2GM 1 /r 
2GM 1 



(d r ®) 2 

(B26) 



and 



Sin = 2tt 



dt I drr 2 

R2 



(dt®) 2 



1 - 2GM to t/r 
2GM to t 



1 



(d r ®) 2 
(B27) 



Here we will make some comments on (|B25|) - (|B27|) . 
First, note that in the limit M\ — >• the upper limit 
of integrate of (|B25p goes to R 2l while the lower limit 
of (IB261) also goes to R 2 . Thus (lB25j) reduces to (44) 
in [27j, while (|B26|) goes to zero. Thus we recover the 
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scenario investigated in [27j. Second, note that in the 
limit M 2 —> 0, B\ = B 2 = B, so that again (|B25|) reduces 
to (44) in [27], while again (|B26j) goes to zero and (|B27D 
reduces to (45) in [27j. Again we recover the scenario 
investigated in [27j . 

We are now in a position to determine the dominate 
parts of the action for the scalar field. First we consider 
the gradient term. We see that in the R 2 — » Rs region 
both the gradient terms in Si and Sjj go to zero since 

— ^ 0, however, the gradient term in Sm goes to a finite 
number, thus we keep the gradient term from (|B27[) . For 
the kinetic term, more care is needed to determine which 
term dominates. For Sm we see that the kinetic term 
is logarithmically divergent, while both Si and Sn are 
(R 2 — Rs)~ X divergent, due to the factor of B. However, 
we also see that the kinetic term in Si goes as (R\ — 
while S11 is logarithmically dependent. Hence, 

since R\ can get arbitrarily close to R^ we can conclude 
that the kinetic term in Si is the dominant term. Hence 
the total action can be written as 



I r B 2 r"'" 



T 



dr ] 



drr 2 (d t $) 2 
i _2GMto 1 



(B28) 



As in [27] one can expand $ in modes and define ma- 
trices A and C that are independent of Ri(t) and R 2 (t), 
here defined by 



C k k' = 47T 



= 4tt / drr 2 f k (r)f k ,(r), (B29) 
Jo 



/ drr 2 I 

Jr s V 



Once again using the standard quantization procedure 
and performing a principle axis transformation, we find 
that for a single eigenmode, the Schrodinger equation 
takes the form 



2m db 2 



*M=i^ (BSD 



where m denotes the eigenvalue of A, b is the eigenmode, 



o 2 (v) = 



K B 2 _ QlB 2 



fh BBi BB\ 
where K denotes the eigenvalue of C and 



7] 



(B32) 



(B33) 



This has identical structure to that found for a single 
massive shell collapsing upon itself and a massive shell 
collapsing onto a PBH. 



Appendix C: Notes on the breakdown of the WKB 
expansion for the Wheeler-de Witt equation 

One implication of our results to which is worth pay- 
ing special attention is related to a famous issue arising 
in first quantization of a scalar field on a collapsing space- 
time [29|, liH EtJ • When one performs such quantization, 
one looks for a wave functional ^ being a solution of the 
Wheeler-de Witt equation ifwdw^ = for the scalar 
field (+ moving shell + gravitational degrees of freedom) 
with vacuum-like initial conditions. The square of this 
wave functional |\I/| 2 then provides the probability den- 
sity to find the shell and emitted Hawking radiation in a 
given state. 

In (3 + l)-dimensions one can look for the solution of 
the Wheeler-de Witt equation in the form ^ = e %s with 
the "action" expanded as 



S = G- 1 S + S 1 + GS 2 - 



(CI) 



where G is the gravitational constant. The expansion 
(|C1|) is nothing but a WKB expansion for the wave func- 
tional. The zeroth order term denoted So is the solu- 
tion of the Einstein-Hamilton- Jacobi equation for the 
background spacetime, the Si term describes effects of 
quantum theory of the scalar field / on the fixed back- 
ground spacetime (it therefore includes the Hawking ef- 
fect) and the wave functional e %Sl satisfies the functional 
Schrodinger equation used in the Sec. El the S 2 term de- 
scribes the backreaction of the background to the Hawk- 
ing radiation in the leading approximation in 1/Mp, etc. 
In the (l + l)-dimensional CGHS case, a similar WKB ex- 
pansion of the wave functional exists: it is either the ex- 
pansion in terms of the dimensionless gravitational con- 
stant or, equivalently, the expansion in terms of the pa- 



f k {r)f k ,{r). (B30) rameter k 
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where TV is the (large) number of 



matter fields / [47j . 

The expansion (|C1|) has been proved to break down 
in the near horizon regime [46|, |47| . This effect seems to 
be puzzling as it cannot be seen in the second-quantized 
version of the theory. Indeed, the standard intuition sug- 
gests that, once the shell is close to its apparent horizon, 
the applicability of the geometric-optics approximation 
for the modes of / improves, and the phases of the modes 
pile up. Correspondingly, one expects that the phase of 
the wave functional \£ oscillates rapidly, i.e, the WKB 
approximation should be better, not worse. 

It is possible to show using the formalism of the quan- 
tum kinetic equation that the WKB breakdown effect 
noticed in [46| can be mostly attributed to the Bose- 
Einstein condensation of "in" excitations and the forma- 
tion of the out-vacuum. As we have discussed, the num- 
ber of instantaneous excitations at late times t — >• + 00 
is given by the solution of the quantum kinetic equation 
([108)) and can be represented in the form of an expansion 
in powers of exp (—t/2Rs): 



N(t, uJk) = N(uj k ) + n(u k )e 



-t/2R s 



+ 0{e- t ' Ra ). (C2) 
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In this expansion, the first term contributes to the polar- 
ization of "in" -vacuum (i.e, formation of the out- vacuum) 
and, correspondingly, to the redefinition of So in (jClj) ). 
The second term proportional to n(uOk) describes the 
Hawking flux (since it corresponds to the finite energy 
flux at spatial infinity, see (jl35ft and subsequent discus- 
sion) and therefore contributes to S\ in (|C1|) . Once the 
vacuum is properly redefined, the zeroth order contri- 



bution ~ Sq in (|C1[) becomes large compared to the first 
order contribution ~ Si , because the second term in (|C2|) 
is subdominant with respect to the first term as t — >> +00. 

We finally note that the WKB breakdown effect cannot 
be entirely attributed to the vacuum polarization, as in 
the slicing we use one expects the quantum gravitational 
effects of "horizon trembling" to become important in 
the near horizon regime [48l-l50( . 



[1] S. W. Hawking, Commun. Math. Phys. 43, 199-220 
(1975) [Erratum-ibid. 46, 206 (1976)]. 
J. D. Bekenstein, Phys. Rev. D7, 2333-2346 (1973). 
S. W. Haw king, Phys. Rev. D14, 2460-2473 (1976). 
J. Preskill, [hep-th/9209058] . 

C. G. Callan, Jr., S. B. Giddings, J. A. Harvey, 
A. Strom inger, Phys. Rev. D45, 1005-1009 (1992). 
[hep-th/91110561 . 

S. M. Christensen and S. A. Fulling, Phys. Rev. D 15, 
2088 (1977). 

T. M. Fiola, J. Preskill, A. Stromi nger, S. P. Triv edi, 
Phys. Rev. D50, 3987-4014 (1994). [hep-th/94031371 . 
A. Ashtekar, F. Preterms, F. M. Ramazanoglu, Phys. 
Rev. D83, 044040 (2011). [arXiv: 1012.00771 [gr-qc]]. 
S. D. Mathur, C lass. Quant. Grav. 26, 224001 (2009). 
[arXiv:0909. 10381 [hep-th]] . 

G. ; t Hooft, Class. Quant. Grav. 16, 3263 (1999) 
|gr-qc/9903084| . 

Y. Aharonov, A. Casher, S. Nussinov, Phys. Lett. B191, 
51 (1987). 

S. B. Giddings, [hep-th/94 12138] , 

S. B. Giddings, Class. Quant. Grav. 28, 025002 (2011). 
[arXiv:0911.3395l [hep-th]]. 
S. B. Giddings. larXiv: 1201. 10371 [hep-th]. 
G. 't Hooft, Nucl. Phys. B 256, 727 (1985). 
L. Susskind, L. Thorlacius and J. Uglum, Phys. Rev. D 
48, 3743 (1993) [hep-th/93060691 . 
A. Strominger, [nep-th/9501071| . 

W. Israel, Z. Yu n, Phys. Rev. D82, 124036 (2010). 
[arXiv:1009.0"879l [hep-th]] . 

S. B. G iddings and W. M . Nelson, Phys. Rev. D 46, 2486 
(1992) [hep-th/9204072] . 

D. Grumiller, W. Kummer and D. V. Vassilevich, Phys. 
Rept. 369, 327 (2002) [hep-th/0204253] , 
L. Susskind, J. Lindesay, "An introduction to black holes, 
information and the string theory revolution: The holo- 
graphic universe," Hackensack, USA: World Scientific 
(2005) 183 p. 

R. Bal binot and A. Fa bbri, Phys. Rev. D 59, 044031 
(1999) [hep-th/9807i23] . 

R. Balbinot, A. Fabbri, V. P. Frolov, P. Nicolini, P. Sut- 
ton and A. Zelnikov, Phys. Rev. D 63, 084029 (2001) 

[hep-th/0012048]. 

P. Townsend, |gr-qc/9707012| 

C. W. Misner, K. S. Thorne and J. A. Wheeler, "Gravi- 
tation," San Francisco 1973, 1279 p. 
[26] J. Ipser, P. Sikivie, Phys. Rev. D30, 712 (1984). 



[27] T. Vachaspati, D. Stoj kovic, L. M. Krauss, Phys. Rev. 

D76, 024005 (2007) gr-qc/0609024] . 
[28] S. A. Fulling, Gen. Rel. Grav. 10, 807 (1979). 
[29] E. Greenwood, D. I. Podolsky and G. D. Starkman, 

JCAP 1111, 024 (2011) [arXiv:1011. 22191 [gr-qc]]. 
[30] P. Carruthers, M. M. Nieto, Rev. Mod. Phys. 40, 411-440 

(1968). 

[31] B. Allen, Phys. Rev. D 32, 3136 (1985). 

[32] E. Mottola, Phys. Rev. D 31, 754 (1985). 

[33] Y. .B. Zeldovich and A. A. Starobinsky, Sov. Phys. JETP 

34, 1159 (1972) [Zh. Eksp. Teor. Fiz. 61, 2161 (1971)]. 
[34] S. Habib, C. Molina-Paris and E. Mottola, Phys. Rev. D 

61, 024010 (2000) [gr-qc/9906120] . 
[35] S. M. Schmidt, D. Blaschke, G. Ropke, S. A. Smolyansky, 

A. V. Prozorke vich, V. D. Tonee v, Int. J. Mod. Phys. E7, 

709-722 (1998) [hep-ph/98092271 - 
[36] Yu. Kagan, "Kinetics of Bose-Einstein condensate forma- 
tion in an interacting Bose gas", in "Bose-Einstein con- 
densation", eds. A. Griffin, D. W. Snoke and S. Stringari, 

Cambridge University Press, 1995. 
[37] N. Kaloper, M. Kleban, A. Lawrence, S. Shenker and 

L. Susskind, JHEP 0211, 037 (2002) [hep-th/020923l] . 
[38] R. Bousso, A. Maloney and A. Strominger, Phys. Rev. D 

65, 104039 (2002) [hep-th/0112218] . 
[39] N. D. Birrell and P. C. W. Davies, Cambridge, Uk: Univ. 

Pr. ( 1982) 340 p. 
[40] L.D. Landau and E.M. Lifshitz, "Mechanics", 

Butterworth-Heinemann (1976) 224 p. 
[41] P. M. Chesler and L. G. Yaffe , Phys. Rev. Lett. 102, 

211601 (2009) [arXiv:0812.2"053l [hep-th]]. 
[42] P. M. Chesler and L. G. Yaffe , Phys. Rev. Lett. 106, 

021601 (2011) [arXiv:1011.3562l [hep-th]]. 
[43] S. Caron-Huot, P. M. Chesler and D. T eaney, Phys. Rev. 

D 84, 026012 (2011) arXiv: 1102. 10731 [hep-th]] . 
[44] P. M. Chesler and D. Teanev. larXiv: 11 12.61961 [hep-th]. 
[45] R. Mic ha and I. I. Tkac hev, Phys. Rev. D 70, 043538 

(2004) [hep-ph/0403"T0T] . 
[46] C. Kiefer, T. P. Singh, Phys. Rev. D44, 1067-1076 

(1991). 

[47] J. -G. Demers, C. Kie fer, Phys. Rev. D53, 7050-7061 

(1996). hep-th/95111471 . 
[48] J. A. Wheeler, "Geometrodynamics," (Academic Press, 

New York, 1962). 
[49] E. Harrison, Phys. Rev. Dl, 2726 (1970). 
[50] Y. W. York, Phys. Rev. D28, 2929 (1983). 



